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Abstract 

There has been significant progress in our understanding of finite-temperature field 
theory over the past decade. In this paper, we review the progress in perturba- 
tive thermal field theory focusing on thermodynamic quantities. We first discuss 
the breakdown of naive perturbation theory at finite temperature and the need for 
an effective expansion that resums an infinite class of diagrams in the perturba- 
tive expansion. This effective expansion which is due to Braaten and Pisarski, can 
be used to systematically calculate various static and dynamical quantities as a 
weak-coupling expansion in powers of g. However, it turns that the weak-coupling 
expansion for thermodynamic quantities are useless unless the coupling constant is 
very small. We critically discuss various ways of reorganizing the perturbative series 
for thermal field theories in order to improve its convergence. These include screened 
perturbation theory (SPT), hard-thermal- loop perturbation theory (HTLPT), the 
^-derivable approach, dimensionally reduced (DR) SPT, and the DR ^-derivable 
approach. 
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1 Introduction 



A detailed understanding of the properties of the deconfined phase of 
QCD is important in several areas of physics. Examples are the evolution of 
the early universe [1], heavy- ion collisions at RHIC (Brookhaven) and LHC 
(CERN) [2,3,4,5,6], and astrophysical objects such as quark stars and super- 
novae [7,8,9]. In both cases, since the QCD running coupling decreases with 
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increasing momentum transfer, the naive expectation is that weak-coupling 
techniques will be sufficient to describe the relevant physics. In the case of 
the RHIC and LHC heavy-ion experiments the expectation is that during the 
collisions that temperatures on the order of 0.2-1 GeV will be generated corre- 
sponding to a strong coupling constant a s ~ 0.2 or g s ~ 2. For cold astrophys- 
ical objects the relevant scale is set by the density which, in the case of quark 
stars, is related to the quark chemical potential. At densities expected at the 
center of neutron stars the quark chemical potential is again on the order of 
0.2-1 GeV and the corresponding strong coupling constant is again a s ~ 0.2. 
Naively one might expect that since a s ~ 0.2 that weak coupling expansions 
could provide a systematic method for calculating observables; however, in 
practice we find that weak-coupling approximations to even thermodynamic 
observables have a small (perhaps vanishing) radius of convergence. This has 
motivated consideration of methods for systematically reorganizing perturba- 
tive expansions in order to obtain more convergent successive approximations 
to observables. In this paper we will review theoretical progress in the per- 
turbative thermodynamics of scalar and gauge theories and provide a critical 
comparison of various methods for improving the convergence of the result- 
ing approximants via non-perturbative variational and effective field theory 
methods. 

In the past 15 years we have improved our understanding of thermal field 
theories considerably. A major problem that spurred progress in the late 1980s 
was the apparent gauge dependence of the gluon damping rate 7. The gluon 
damping rate had been calculated in ordinary perturbation theory by a number 
of authors in various gauges [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]. 
It turned out that the result obtained for the damping rate depended on the 
gauge-fixing condition. Even the sign of the damping rate seemed to depend 
on the gauge which would be worrisome since, in those gauges where the sign 
is negative, this indicated a plasma instability. However, since the damping 
rate is a physical quantity, being the inverse of the lifetime of a collective 
excitation in the plasma, it cannot be gauge dependent. The problem was 
solved by Pisarski [26], who pointed out that the one-loop calculations which 
had been performed to date were incomplete and that in order to obtain the 
correct leading-order result in the coupling g, one must sum an infinite subset 
of diagrams. A detailed account of the problem has been given in Ref . [27] . The 
general program of resummation was later developed by Braaten and Pisarski 
in Ref. [28], where they showed that effective propagators and vertices must 
be used whenever the momentum characterizing the external lines are of order 

gT. 

The next significant step was made a few years later, when perturbative 
calculations of the free energy for high-temperature field theories were pushed 
to higher order in the coupling. The free energy to order g 4 has been calcu- 
lated in massless 4 -theory in Ref. [29], for QED in Ref. [30], and for QCD in 
Ref. [31]. In particular the calculations by Arnold and Zhai [31] were impor- 
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tant since they were completely analytic. Since g 4 is the first order at which 
renormalization effects enter, these corrections are important beyond deter- 
mining another term in the perturbative expansion. Calculations were later 
pushed to order g 5 in Ref. [32,33] for $ 4 -theory, in Refs. [34,35] for QED, and 
in Refs. [36,37] for nonabelian gauge theories. 

If one is interested in static quantities such as the pressure or screen- 
ing lengths, there exists an efficient alternative to explicit resummation. It 
is based on the idea of dimensional reduction [38,39,40,41] and effective field 
theory methods [43]. The basic idea is that equilibrium properties of a 3+1 di- 
mensional field theory at high temperature can be calculated using an effective 
field theory in three spatial dimensions. This follows from the fact that in the 
imaginary time formalism, the Matsubara frequencies act like masses in the 
propagators and the nonstatic modes decouple according to the Appelquist- 
Carrazone theorem [44]. This approach was developed into tool for quanti- 
tative calculations by Braaten and Nieto [33,37], and by Farakos, Kajantie, 
Rummukainen, and Shaposhnikov [45]. 

However, these results showed that the weak-coupling expansion for the 
pressure of a thermal field theory only converges if the coupling constant is 
very small. As successive terms in the weak-coupling expansion are added, 
the predictions change wildly and the sensitivity to the renormalization scale 
/i grows. In the case of QCD, the order-g 3 contribution is smaller than the 
order-g 2 contribution only if a s < 1/20, which corresponds to a temperature 
of 10 5 GeV. This is orders of magnitude higher than the temperature expected 
to be generated at RHIC and LHC. It is clear that a reorganization of the 
perturbation series is essential if perturbative calculations are to be of any 
quantitative use at temperatures accessible in heavy-ion collisions. 

There are some proposals for reorganizing perturbation theory in hot field 
theories that are basically mathematical manipulations of the weak coupling 
expansion. The methods include Pade approximates [46], Borel resummation 
[47] , and self-similar approximates [48] . Additionally we note that in Ref. [49] 
Pade and Borel resummation techniques have been applied separately to the 
soft and hard contributions to the QCD free energy. The methods listed above 
can be used to construct more stable sequences of successive approximations 
which agree with the weak-coupling expansion when expanded in powers of 
g. However, there are two major limitations of this approach: the first is that 
these methods can only be applied to quantities for which several orders in 
the weak-coupling expansion are known, so they are limited in practice to the 
thermodynamic functions; the second is that selecting which Pade approxi- 
mant is used is a non-trivial task with different Pade approximants sometimes 
giving dramatically different results with the difference between resulting ap- 
proximants being on the order of the variations of the successive perturbative 
approximations. 

The free energy can also be calculated nonperturbatively using lattice 
gauge theory [50]. The thermodynamic functions for pure-glue QCD have been 
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calculated with high precision by Boyd et al [51]. There have also been cal- 
culations with Nf = 2 and 4 flavors of dynamical quarks [52]. These lattice 
studies find that the free energy is very close to zero near T c and, as the 
temperature increases, the free energy increases rapidly from T c to 2T C and 
then slowly approaches that of an ideal gas of massless quarks and gluons. 
Continued advances in lattice gauge theories will provide even more precise 
information about the QCD phase transition. In practice, however, lattice 
gauge theories have some limitations. For instance, the method can only be 
applied to static quantities and many of the more promising signatures for a 
quark-gluon plasma involve non-equilibrium dynamical quantities. One would 
therefore like a coherent framework for calculating both static and dynamical 
quantities that would reproduce lattice data. A second major limitation of lat- 
tice methods is that the Monte Carlo approach used fails at nonzero baryon 
number density due to the sign problem of the action. It is therefore impor- 
tant to have perturbative and non-perturbative analytical techniques which 
can be used in this region of the QCD phase diagram. Note, however, that 
recent advances have made it possible to use lattice methods to calculate the 
pressure at small baryon densities [53,54,55,56,57,58,59,60,61]. 

One analytical approach to the thermodynamics of QCD is the application 
of so-called quasi-particle models [62,63,64]. In these models, the quark-gluon 
plasma is simply an ideal gas of quasi-particles with temperature-dependent 
masses. The zero-point energy of the quasi-particles is replaced by a bag con- 
stant which is necessary to enforce thermodynamic consistency. By fitting the 
parameters to finite-temperature lattice data one can obtain a quite reasonable 
description of the thermodynamics of QCD. In addition, using thermodynamic 
self-consistency results obtained at T = can be extended to nonzero chem- 
ical potential [65] (see also Ref. [66]). While these approaches seem to do a 
good job at fitting the existing lattice data, it would be preferable to have a 
framework which would incorporate the same physics, namely quasi-particle 
degrees of freedom, in a more systematic way. 

One way to accomplish this goal is to reorganize the weak-coupling ex- 
pansion based on a variational approach. The free energy T is then expressed 
as the variational minimum of a thermodynamic potential Q(T, a s ; m 2 ) that 
depends on one or more variational parameters that are collectively denoted 
by m 2 . One such approach is the ^-derivable approximation in which the ex- 
act propagator is used as an infinite set of variational parameters [67,68,69]. 
The variational principle leads to a gap equation for the self-energy which 
is, in general, extremely difficult to solve. Only in simple cases where the 
self-energy is independent of the external momentum, is the theory easy to 
solve. Approximate solutions of the two-loop ^-derivable approximation in 
terms of HTL self-energies have been obtained for QCD at finite tempera- 
ture [70,71,72,73] and chemical potential [74,75]. Additionally, there has been 
progress recently both with respect to important issues such as renormal- 
ization [76,77] and gauge (in-) dependence [78,79] as well as the development 
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of calculational technology [80]. The ^-derivable approach has recently been 
reviewed in Refs. [81] and [82]. 

The intractability of the ^-derivable approach motivates the use of sim- 
pler variational approximations. One such strategy that involves a single vari- 
ational mass parameter m has been called optimized perturbation theory [83], 
variational perturbation theory [84], or the linear 5 expansion [85]. This strat- 
egy was applied to the thermodynamics of the massless 4 field theory by 
Karsch, Patkos and Petreczky under the name screened perturbation theory 
[86]. The method has also been applied to spontaneously broken field the- 
ories at finite temperature [87]. The calculations of the thermodynamics of 
the massless 4 field theory using screened perturbation theory have been ex- 
tended to three loops [88]. The calculations can be greatly simplified by using 
a double expansion in powers of the coupling constant and m/T [89]. Screened 
perturbation theory has also been generalized to gauge theories [90,91] and 
is then called hard-thermal-loop (HTL) perturbation theory. Two-loop calcu- 
lations in pure-glue QCD [92] and QCD with dynamical fermions have been 
performed [93]. 

The dimensionally reduced effective field theory can also be used as the 
starting point for nonperturbative calculations of static properties. The coef- 
ficients in the effective Lagrangian are calculated using perturbation theory, 
but calculations in three dimensions are performed nonperturbatively using 
analytic methods and/or lattice gauge theory. Unfortunately, dimensional re- 
duction has similar limitations to ordinary lattice gauge theory: it can be 
applied only to static quantities and only for baryon number densities which 
are much smaller than the temperature. Unlike ordinary lattice gauge the- 
ory, however, light dynamical quarks do not require any additional computer 
power, because they only enter through the perturbatively calculated param- 
eters of £ e fr- This method has been applied to the Debye screening mass for 
QCD [94] as well as the pressure [95,96]. 

The paper is organized as follows. In Sec. 2, we discuss the need for re- 
summation and the Braaten-Pisarski resummation program. In Sec. 3, we dis- 
cuss the weak-coupling expansion of the thermodynamic quantities in detail. 
Sec. 4 is devoted to dimensional reduction and the application of effective field 
theory methods to hot field theories In Sees. 5 and 6, we discuss screened per- 
turbation theory and hard-thermal-loop perturbation theory. The ^-derivable 
approach is reviewed in Sec. 7. In Sees. 8 and 9 we cover dimensionally re- 
duced screened perturbation theory and dimensionally reduced ^-derivable 
approaches. In Sec. 10, we summarize and conclude. There are three appen- 
dices where we list sum-integrals and integrals needed as well as some technical 
details of the calculations. 
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2 The Need for Resummation 



In this section, and in the rest of the paper, we consider thermal field 
theories at high temperatures, which means temperatures much higher than 
all zero-temperature masses or any mass scales generated at zero temperature. 

It has been known for many years that naive perturbation theory, or the 
loop expansion breaks, down at high temperature due to infrared divergences. 
Diagrams which are nominally of higher order in the coupling constant con- 
tribute to leading order in g. A consistent perturbative expansion requires the 
resummation of an infinite subset of diagrams from all orders of perturbation 
theory. We discuss these issues next. 



2. 1 Scalar field theory 



We start our discussion by considering the simplest interacting thermal 
field theory, namely that of a single massless scalar field with a -interaction. 
The Euclidean Lagrangian is 



Perturbative calculations at zero temperature proceed by dividing the La- 
grangian into a free part and an interacting part according to 

£free=i(^) 2 , (2) 

A„t = f^> 4 . (3) 

Radiative corrections are then calculated in a loop expansion which is equiv- 
alent to a power series in g 2 . We shall see that the perturbative expansion 
breaks down at finite temperature and the weak-coupling expansion becomes 
an expansion in g rather than g 2 . 

We will first calculate the self-energy by evaluating the relevant diagrams. 
The Feynman diagrams that contribute to the self-energy up to two loops are 
shown in Fig. 1. 



1-loop : 2-loop 



Fig. 1. One- and two-loop scalar self-energy graphs. 

The one-loop diagram is independent of the external momentum and the 
resulting integral expression is 



P 2 



6 



9 2 rp2 

24 



: m 



(4) 



where the superscript indicates the number of loops. The sum-integral Jl p 
represents a summation over Matsubara frequencies and integration of spatial 
momenta in d = 3— 2e dimensions. It is defined in Eq. (A.l) * . The above sum- 
integral has ultraviolet power divergences that are set to zero in dimensional 
regularization. We are then left with the finite result (4), which shows that 
thermal fluctuations generate a mass for the scalar field of order gT. This 
thermal mass is analogous to the Debye mass which is well-known from the 
nonrelativistic QED plasma. 

We next focus on the two- loop diagrams and first consider the double- 
bubble in Fig. 2 (b). 

_Q_= XI + _8_ + i£ + - 

(a) (b) (c) 

Fig. 2. Bubble diagrams contributing to the scalar self-energy. 

This diagram is also independent of the external momentum and gives the 
following sum-integral 



This integral is infrared divergent. The problem stems from the loop with two 
propagators. In order to isolate the source of the divergence, we look at the 
contribution from the zeroth Matsubara mode to the Q integration 

q 

which is quadratically infrared divergent. This infrared divergence indicates 
that naive perturbation theory breaks down at finite temperature. However, 
in practice this infrared divergence is screened by the thermally generated 
mass and we must somehow take this into account. The thermal mass can be 
incorporated by using an effective propagator: 

AK,p)= - 1 - , (7) 



P 2 + m 



with m ~ gT <C T. 



* For an introduction to thermal field theory and the imaginary time formalism see 
Refs. [97] and [98]. 
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If the momenta of the propagator is of order T or hard, clearly the thermal 
mass is a perturbation and can be omitted. However, if the momenta of the 
propagator is of order gT or soft, the thermal mass is as large as the bare in- 
verse propagator and cannot be omitted. The mass term in the propagator (7) 
provides an infrared cutoff of order gT. The contribution from (6) would then 
be 

" \4rh T I W^f- = (9 GL) + ° ^ ' <8) 

Since m ~ gT, this shows that the double-bubble contributes at order g 3 T 2 to 
the self-energy and not at order g 4 T 2 as one might have expected. Similarly, 
one can show that the diagrams with any number of bubbles like Fig. 2c are 
all of order g 3 . Clearly, naive perturbation theory breaks down since the order- 
g 3 correction to the thermal mass receives contributions from all loop orders. 
On the other hand, the three-loop diagram shown in Fig. 3, is of order g A T 2 
and thus subleading. Therefore, we only need to resum a subset of all possible 
Feynman graphs in order to obtain a consistent expansion in g. 



Fig. 3. Subleading three-loop self-energy diagram. 



If we use the effective propagator to recalculate the one-loop self-energy, 
we obtain 



2 9 



T 



+ m? 
1 

p2 _|_ m 2 



1 



,p2 



O (m 



ff T 2 

24 



4tt 



+ O (g 2 ) 



(9) 



where here, and in the following, the prime on the sum-integral indicates that 
we have excluded the n = mode from the sum over the Matsubara frequen- 
cies. The order g 3 corresponds to the summation of the bubble diagrams in 
Fig. 2, which can be verified by expanding the effective propagator (7) around 
m — 0. Thus by taking the thermal mass into account, one is resumming an 
infinite set of diagrams from all orders of perturbation theory. 

The self-energy (4) is the first example of a hard thermal loop (HTL). Hard 
thermal loops are loop corrections which are g 2 T 2 /P 2 times the corresponding 
tree-level amplitude, where P is a momentum that characterizes the external 
lines. From this definition, we see that, whenever P is hard, the loop correction 
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Fig. 4. One-loop photon self-energy diagram. 



is suppressed by g 2 and is thus a perturbative correction. However, for soft 
P, the hard thermal loop is 0(1) and is therefore as important as the tree- 
level contribution to the amplitude. These loop corrections are called "hard" 
because the relevant integrals are dominated by momenta of order T. Also 
note that the hard thermal loop in the two-point function is finite since it is 
exclusively due to thermal fluctuations. Quantum fluctuations do not enter. 
Both properties are shared by all hard thermal loops. 

What about higher-order n-point functions in scalar thermal field theory? 
One can show that within scalar theory the one-loop correction to the four- 
point function for high temperature behaves as [27] 

r( 4 Wiog(T/ P ) , (io) 

where p is the external momentum. Thus the loop correction to the four-point 
function increases logarithmically with temperature. It is therefore always 
down by one power of g, and it suffices to use a bare vertex. More gener- 
ally, it can be shown that the only hard thermal loop in scalar field theory is 
the tadpole diagram in Fig. 1 and resummation is taken care of by including 
the thermal mass in the propagator. In gauge theories, the situation is much 
more complicated as we shall see in the next section. 

2.2 Gauge theories 

In the previous section, we demonstrated the need for resummation in 
a hot scalar theory. For scalar theories, resummation simply amounts to in- 
cluding the thermal mass in the propagator and since the running coupling 
depends logarithmically on the temperature, corrections to the bare vertex are 
always down by powers of g 2 . In gauge theories, the situation is more com- 
plicated. The equivalent HTL self-energies are no longer local, but depend 
in a nontrivial way on the external momentum. In addition, it is also nec- 
essary to use effective vertices that also depend on the external momentum. 
It turns out that all hard thermal loops are gauge- fixing independent. This 
was shown explicitly in covariant gauges, Coulomb gauges, and axial gauges. 
They also satisfy tree-level like Ward identities. Furthermore, there exists a 
gauge invariant effective Lagrangian, found independently by Braaten and Pis- 
arski [99] and by Taylor and Wong [100], that generates all of the hard thermal 
loop n-point functions. From a renormalization group point of view this is an 
effective Lagrangian for the soft scale gT that is obtained by integrating out 
the hard scale T. We return to the HTL Lagrangian in Section 6. 
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2.2.1 Polarization tensor 

We next discuss in some detail the hard thermal loop for the polarization 
tensor H^. For simplicity, we discuss QED. The Feynman diagram for the one- 
loop self-energy is shown in Fig. 4 and results in the following sum-integral 



IV(w„,p) = e 2 i Tr 

T{K} 



K 2 (P + K) 2 



where Tr denotes the trace over Dirac indices. After taking the trace, the 
self-energy becomes 



1 

{K}JO 



+25 llu P 2 e 2 rf; 



1 



J{K} K 2 [P + K) 2 



P^ + PyKp 



{K} K 2 (P + Kf 



,(12) 



where we have assumed, for now, that d = 3. 

We first consider the spatial components of n MJ/ (o; n , p). The sum over Mat- 
subara frequencies can be written as a contour integral in the complex energy 
plane. After analytic continuation, we obtain 



ILij(w,p) = e 2 



duo 
2tH 



Akjkj - (uj 2 - p 2 )5jj + 2(pikj + kjpj) 
(k 2 -uj 2 )[(p + k) 2 -(uj + uj ) 2 } 



25i 



(k 2 - uj 2 ) 



tanh 



(3uj 



(13) 



In order to extract the hard thermal loop, we notice that terms that contain 
one or more power of the external momentum can be omitted since the external 
momentum p is assumed to be soft. The self-energy then becomes 



(k 2 -uj 2 )[(p + k) 2 -(uj + uj ) 2 } 



tanh 



(3uj 



(k 2 - uj 2 ) 

After integrating over the energy uj , we obtain 

lUji^p) = -2e% J i (1 - 2n F (k)) + 2e 2 J 
k k 

x j (1 -n F (k) -n F (|P + k D) 
— [np(k) — np(\p + k|)] 



(14) 



k\p + k| 

+ 



1 



k + |p + k| + uj k + |p + k| — uj 
1 1 



+ 



p + k| — k + uj |p + k| — k — uj 



,(15) 
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where uf(x) = l/(exp(/3x) + 1) is the Fermi-Dirac distribution function. The 
zero-temperature part of Eq. (15) is logarithmically divergent in the ultravio- 
let. This term depends on the external momentum and is cancelled by standard 
zero-temperature wavefunction renormalization. We next consider the terms 
that depend on temperature. In the case that the loop momentum is soft, the 
Fermi-Dirac distribution functions can be approximated by a constant. The 
contribution from the integral over the magnitude of k is then of order g 3 and 
subleading. When the loop momentum is hard, one can expand the terms in 
the integrand in powers of the external momentum. We can then make the 
following approximations 

n F (\p + k|) «n F (fc) + ^^P-k , (16) 

(17) 



|p + k| + k±u 2k' 



|p + k| - k±u p-k±cu ' 

where k = k//c is a unit vector. Thus the angular integration decouples from 
the integral over the magnitude k. This implies 

. . 4e 2 7 ,, , /, \ [ d£l cu 

Il ij (u,p) = — / dkkn F (k) / -kikj , 

7T J J An uj - p-k 

2 r d£l uj f f 
= m D / -: fhkj , (19) 

where we have defined m 2 D = e 2 T 2 /3. 

The other components of the self-energy tensor Il M ^(cJ,p) are derived in 
the same manner. One finds [98] 

n^.p)^/^-^-/) , (20) 

U ,(uJ,p)=m 2 D (21) 

In d dimensions, we can compactly write the self-energy tensor as 

W v {p) = m 2 D [T^{p, -p) - n»n v \ , (22) 

where n specifies the thermal rest frame is canonically given by n — (1, 0), we 
have defined 



m 2 



^-4( d -l)e^ ±, (23) 
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and the tensor T^ u (p, q), which is defined only for momenta that satisfy p+q = 
0, is 



(24) 



The angular brackets indicate averaging over the spatial directions of the light- 
like vector y = (l,y). The tensor T^ u is symmetric in // and v and the self- 
energy (22) satisfies the Ward identity: 



(25) 



Because of this Ward identity and the rotational symmetry around the p- 
axis, one can express the self-energy in terms of two independent functions, 
II T (u;,p) and U L (u,p): 



IV (w, p) = n L (w, p) f — 



+ 



cj 2 — p 2 
U t (uj, p) — U L (u, p) 



(Sij - pipj) g ju , (26) 



where the functions n T (u;,p) and U L (u : p) are 



U T (uj,p) = -(Sij -PiPj)Uij(u,p) , (27) 
U L (co,p) = -Uoo(uJ,p) . (28) 

In three dimensions, the self-energies Hr(c<;,p) and Ul(uj,p) reduce to 



U T (uj,p) 



m 2 D uj 2 



2 p 2 
U L (uj,p)=m 2 D 



p 2 -uj 2 uj + p 

1 + — log 

loop uj — p 



UJ UJ+p 

1 - — log 

2p uj — p 



(29) 
(30) 



The hard thermal loop in the photon propagator was first calculated by Silin 
more than forty years ago [101]. The hard thermal loop in the gluon self-energy 
was first calculated by Klimov and Weldon [102,103]. It has the same form as 
in QED, but where the Debye mass is replaced by 



2 2 

m D = g 



(31) 



where Ca = N c is the number of colors and Nf is the number of flavors. When 
d = 3 the QCD gluon Debye mass becomes 



ml = 1(Ca + \N f ) g 2 T 2 



(32) 
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2.2.2 Fermionic self- energy 

The electron self-energy is given by 

E(P) = m) lti r^) , (33) 
where 

T"(p) = /-^\ , (34) 

and m/ is the thermal electron mass 

m} = -3e 2 f -L. (35) 
In QCD, the quark mass is given by 



"*«=- 3C «2v>^' (36) 

-2. 5 Higher n-point functions 
In gauge theories, there are also hard thermal loops involving vertices. 
For instance, the one-loop correction to the three-point function in QED, can 
compactly be written as 

r^ = y -m 2 f f>*(p,q,r) , (37) 
where the tensor in the HTL correction term is only defined for p — q + r = 0: 

f*(p,q,r) = (y* (—*—)) . (38) 

The quark-gluon vertex satisfies the Ward identity 

p^(p,q,r) = S- 1 (q)-S- 1 (r) } (39) 

where S(q) is the resummed effective fermion propagator. 

In QED there are, in fact, infinitely many amplitudes with hard thermal 
loops. To be precise, there are hard thermal loops in all n-point functions with 
two fermion lines and n — 2 photon lines. In nonabelian gauge theories such 
as QCD, there are in addition hard thermal loops in amplitudes with n gluon 
lines [28]. 

3 The Weak-coupling Expansion 

The Braaten-Pisarski resummation program has been used to calculate the 
thermodynamic functions as a weak-coupling expansion in g. They have now 
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been calculated explicitly through order g 5 for massless c/> 4 theory [32,33], for 
QED [30,34,35], and for nonabelian gauge theories [31,36,37]. In this section, 
we review these calculations in some detail. 



3.1 Scalar field theory 

The simplest way of dealing with the infrared divergences in scalar field 
theory is to reorganize perturbation theory in such a way that it incorpo- 
rates the effects of the thermally generated mass m into the free part of the 
Lagrangian. One possibility is to divide the Lagrangian (1) into free and in- 
teracting parts according to 

C lTec = l -{d,<P) 2 + l -m 2 <p\ (40) 

C mt = Cp'- l -m^\ (41) 

Both terms in Eq. (41) are treated as interaction terms of the same order, 
namely g 2 . However, the resummation implied by the above is rather cumber- 
some when it comes to calculating Green's function with zero external energy. 
Static Green's functions can always be calculated directly in imaginary time 
without having to analytically continue them back to real time. This implies 
that we can use Euclidean propagators with discrete energies when analyz- 
ing infrared divergences which greatly simplifies the treatment. In particular, 
since only propagators with zero Matsubara frequency have no infrared cutoff 
of order T only for these modes is the thermal mass of order gT relevant as an 
IR cutoff. Thus, another possibility is to add and subtract a mass term only 
for the zero-frequency mode. This approach has the advantage that we do not 
need to expand the sum-integrals in powers of m 2 /T 2 in order to obtain the 
contribution from a given term in powers of g 2 . We will follow this path in the 
remainder of this section and write 

£ fr cc = i(<V>) 2 + l -m 2 <p 2 5 P0fl , (42) 
C^= g ^/-\m 2 <P% 0fl . (43) 
The free propagator then takes the form 

+ ( 44 ) 

This way of resumming is referred to as static resummation. It is important 
to point out that this simplified resummation scheme can only be used to cal- 
culate static quantities such as the pressure or screening masses. Calculation 
of dynamical quantities requires the full Braaten-Pisarski resummation pro- 
gram. The problem is that the calculation of correlation functions with zero 
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J~Q J~la. J~ \h J--1& J~2h J~2c J~ '2d 

Fig. 5. Diagrams which contribute up to three-loop order in scalar perturbation 
theory. A boldfaced x indicates an insertion of m 2 . 

external frequencies cannot unambiguously be analytically continued to real 
time [116]. 

We next consider the calculation of the free energy through order g 5 in 
scalar field theory. This involves the evaluation of vacuum graphs up to three- 
loop order. The diagrams are those shown in Fig. 5. 



3.1.1 One Loop 

The one-loop contribution to the free energy is 

Fo = \t J log {f + m 2 ) + ^ log P 2 , 



(45) 



Using the integrals and sum-integrals contained in the appendices the result 
for this diagram in the limit e — > is 



7r 2 . Tm 3 

To = T 4 

90 12vr 



(46) 



3.1.2 Two Loops 

The two-loop contribution to the free-energy is given by 



lb 



(47) 



with 



jF lb = —-m 2 T 



2 J p ' z + rri 1 
v 



(48) 
(49) 



The result for these diagrams in the limit e — > is 



^ vr 2 T 4 [5 



1 + e 4 + 4 



C'(-i)' 



5^ 



-a 



3/2 



1 + e 4 + 2 



C(-i) 
C(-i) ' 

c(-i), 



UttIV 



15 



ib : 



tt 2 T 4 5x/6 
90 2 



-a 



3/2 



(51) 



where we have kept all terms that contribute through order e, because they 
are needed for the counterterm diagrams in the three-loop free energy. 

3.1.3 Three Loops 

The three-loop contribution is given by 



?2 = ^2a + ?2h + ^2c + ^2d + 1 ~^A 1 g Z , 



(52) 



where the expressions for the diagrams are 



l j:' l 

+ 



l 



^2b = -—g A zL 



4' 



p2 _|_ m 2 TpP 2 
111 ] 



T 



(p 2 + m 2 y 



48^ TpqrP 2 Q2R2(P + Q + R)2 

LgW I' - J 1 L 

48 J v 2 + rn 2 q 2 + m 2 r 2 + : 



pgr 



p2 _|_ m 2 g2 _|_ m 2 r 2 _|_ m 2 (p + q + r) 2 + m? 



+ 0(g 6 ),(M) 



^ 2c ~ ( T / (p2 + m?) + ^pP 2 ) ( T / (p 2 + m 2 ) 2 ) ' 



i J (p 2 + m 2 ) 2 

p 



(55) 
(56) 



In the appendix, we show how to rewrite the contribution from jF 2b . The result 
for these diagrams in the limit e — > is 



J~\ 2a — 



2b 



•7~2c : 
^2d : 



^ |_ We /2 _ 5 

90 I 8 8 
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V47r77 
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V47rTy 



6e~ 
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7T 2 T 4 

90 



5^ 



a 3 / 2 - -a 2 



7T 2 T 4 5y^ 
90 
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(57) 



(58) 
(59) 
(60) 
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3.1.4 Pressure through g 5 

Combining the one- , two-, and three-loop contributions given by Eqs. (45), 
(47), and (52), respectively, gives the free energy through order g 5 



2t4 



7C Z T 

~9CJ 



5VQ 



-a 



3/2 



15 / a 59 

lOET 

4 V 2nT 15 



-31og2 + 7 



C(-i 

5 



C(-3) 



15y/E r u 2 , 5 

log — — log a + - 

2 V 2ttT 3 5 6 



2, 1 2 C (-1 

- log2 + - log 3 + -7 - -^7 — { 
3 & 3 & 3' 3C-1 



a 



5/2 



(61) 



The pressure through order g 5 was first calculated using resummation by Par- 
wani and Singh [32] and later by Braaten and Nieto using effective field theory 
[33]. 

The renormalization group equation for the coupling g 2 is 



da , 
a— = da 



(62) 



Using Eq. (62), one can verify that the free energy (61) is RG-invariant up to 
corrections of order g 6 log g. 
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Fig. 6. Weak-coupling expansion for the pressure to orders g 2 , g 3 , g 4 , and g 5 nor- 
malized to that of an ideal e function of g(2irT). 

In Fig. 6, we show the successive perturbative approximations to the pres- 
sure as a function of g(2irT). The bands are obtained by varying the renormal- 
ization scale fi by a factor of two around the central value /i = 2tcT. The lack 
of convergence of the weak-coupling expansion is evident from this Figure. 
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The band obtained by varying // by a factor of two is not necessarily a good 
measure of the error, but it is certainly a lower bound on the theoretical error. 
Another indicator of the theoretical error is the deviation between successive 
approximations. We can infer from Fig. 6 that the error grows rapidly when 
g{2^T) exceeds 1.5. 



3.2 Gauge theories 



In this section, we discuss the application of resummation techniques 
to QCD. The Euclidean Lagrangian for an SU(N C ) gauge theory with Nf 
fermions in the fundamental representation is 

£ = \g; u G% + ^D^, (63) 

where G a = d^A^ — d v A a + gf abc A b A^ is the field strength, g is the gauge 
coupling and f abc are the structure constants. The covariant derivative is = 
<9 M + igA a ^T a , where T a are the generators in the fundamental representation. 
In terms of the group-theory factors Ca, C f , d A , d F , and T F , we have the 
relations 

rabc cabd r~t red 

J J =C A d , 

5 aa = d A , 
(T a T% = C F 5 i3 , 

Tr (T a T b )=T F 5 ab , 



S u = d F . (64) 

Moreover, for SU(N C ), we have C A = N c , C F = (N? - 1)/2N C , T F = N f /2 
d A = Nl - 1, and d F = N c N f . 

If we are only interested in static quantities, we can apply the same sim- 
plified resummation scheme also to gauge theories. Thus we are interested in 
the static limit of the polarization tensor U^ u (uj n , p). In that limit U F van- 
ishes and Il £ = m 2 D . In analogy with the scalar field theory, we rewrite the 
Lagrangian by adding and subtracting a mass term |m|,AQAg(5p 0i o- One of the 
mass terms is then absorbed into the propagator for the timelike component 
of the gauge field Aq, while the other is treated as a perturbation. 

The free energy through g 5 requires the evaluation of diagrams up to three 
loops. The strategy is the same in the scalar case, where one distinguishes 
between hard and soft loop momenta. The result for QCD (iV c = 3) is 

87r 2 T 4 f 21 T 15/ 5 A \a, / 1, T \ 3/2 fa s \ 3 / 2 

T= <1 + —Nf [1 + —Nf)— + 30 [1 + -Nf) — ' 

45 I 32 } 4 V 12 V 7T V 6 V V vr 



+ 



2 135 /. 1 



237.2 + 15.97JV/ - 0.413iV; + — ( 1 + -N f ) log 
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g(27tT) 



Fig. 7. Weak-coupling expansion for the pressure of SU(3) gauge theory with Nf = 2 
to orders g 2 , g 3 , g 4 , and g 5 normalized to that of an ideal gas as a function of g(2nT). 
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The free energy for QCD through order g 4 was first derived by Arnold 
and Zhai [31]. Later it was extended to order g 5 by Zhai and Kastening [36] 
using the above resummation techniques, and by Braaten and Nieto using 
effective field theory [37]. We shall return to that in the next section. The 
order-g 5 contribution is the last contribution that can be calculated using 
perturbation theory. At order g 6 , perturbation theory breaks down due to 
infrared divergences in the magnetic sector [104]. 

In Fig. 7, we show the weak-coupling expansion for the pressure of SU (3) 
gauge theory with Nf = 2 to orders g 2 , g 3 , g 4 , and g 5 divided by that of 
an ideal g function of g(2irT). The fluctuations of the predictions as 

one adds terms in the weak-coupling expansion is evident from the Figure. 
Additionally, if the series were well-behaved, the expectation would be that 
the scale variation would be reduced as each additional order in g was included; 
however, as can be seen from Fig. 7, for low-temperatures, the variation of the 
pressure obtained by varying fi by a factor of two is becoming larger with each 
additional perturbative order. This points to the necessity to resum higher logs 
in the final result. As we will see in Sec. 4 this is possible using dimensional 
reduction and effective field theory methods. 
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4 Dimensional Reduction 



In the imaginary time formalism, bosonic fields are periodic and fermionic 
fields are antiperiodic in the time direction with period f3. A generic field can 
therefore be expanded in a Fourier series: 



where u n = 2nnT for bosons and uj n = (2n + V)irT for fermions. The free 
propagator of the n-th component of the field is given by 



One can therefore view a thermal field theory in d + 1 dimensions as 
a Euclidean field theory in d dimensions with infinitely many fields, where 
the Matsubara frequencies act as masses. The nonzero bosonic modes and 
the fermionic modes have a mass of order T, while the zero-frequency mode is 
massless. If we would use a resummed propagator the propagator for the n = 
Matsubara mode would include an effective mass of order gT, while the masses 
for the nonzero bosonic modes and the fermionic modes would still be of order 
T. If the coupling constant g is small, the scales T and gT are well separated 
and one can use effective field theory methods to organize the calculations 
of physical quantities such that one takes care of one scale at a time. Not 
only can effective field theory methods be used to unravel contributions from 
various momentum scales, they also simplify practical calculations by treating 
one scale at a time. More specifically, one constructs a three-dimensional field 
theory for the n = bosonic mode by integrating out all the nonzero bosonic 
modes as well as all the fermionic modes. This process is called dimensional 
reduction (DR) and the nonzero Matsubara modes decouple [38,39,40,41,42]. 

The simplest way to construct the effective field theory for the zero- 
frequency bosonic modes is to use methods from effective field theory [43]. 
Having identified the degrees of freedom and their symmetries in the effec- 
tive three-dimensional theory, one writes the most general Lagrangian C cS 
for the fields that satisfies those symmetries. The effective Lagrangian has 
infinitely many interaction terms and is, in general, nonrenormalizable. The 
terms in C e g are not arbitrary but restricted by the symmetries of C e g and 
the coupling constants are determined by a matching procedure. One calculates 
static correlators in the full and effective theory and demands that they be 
the same at long distances scales R 3> l/T [33]. Unfortunately, the matching 
procedure is complicated by the occurrence of ultraviolet divergences. There 
are two sources of UV divergences. One source is associated with the original 
four-dimensional field theory and these divergences are removed by standard 
renormalization at zero temperature. The other source is associated with UV 




(66) 



20 



divergences in the effective theory in three dimensions and these are gener- 
ated by integrating out the nonzero Matsubara modes. Thus the parameters 
of the effective theory must depend on the ultraviolet cutoff A that is used to 
regulate the loop integrals in the effective theory in such a way that physical 
quantities do not depend on A. The ultraviolet cutoff A plays the role of a 
factorization scale that separates the momentum scale T from the momentum 
scale gT which can be described by the effective three-dimensional theory. 

Modern renormalization theory implies that the static correlators at dis- 
tances R^> 1/T in the full theory can be reproduced to any desired accuracy 
provided we include sufficiently many operators in £ e fr and tune their coef- 
ficients as functions of the coupling constants in the full theory, T, and A 
[33]. The matching procedure is carried out using strict perturbation theory, 
which is identical to naive perturbation theory of Sec. 2. This expansion is, 
of course, plagued with the infrared divergences due to the long-range forces 
mediated by massless particles. We know these divergences are screened, but 
taking screening into account requires summing an infinite set of higher-order 
diagrams. However, we can still use strict perturbation theory to determine 
the coefficients of C e s because they encode the physics on the scale T and 
are insensitive to the physics on the scale gT. If we make the same incorrect 
assumptions about screening in the effective theory, the infrared divergences 
will be the same in the correlators that are being calculated in the two theories 
and they will cancel in the matching procedure [33]. 

4-1 Scalar field theory 

In this subsection, we discuss dimensional reduction and calculations in the 
effective theory for a thermal scalar field theory defined by the Lagrangian (1). 
In particular, the weak-coupling expansion for the pressure (61) is rederived. 

4-1.1 Effective Lagrangian 

The three-dimensional effective field theory obtained by dimensional re- 
duction describes a scalar field, 0(x), that up to normalization can be identified 
with the zero-frequency mode of the original field: 



Based on the general discussion above, we can easily write down the effective 
Lagrangian for a massless scalar field theory. There is a three-dimensional 




(68) 



o 



rotational symmetry as well as the symmetry 



— 0. One finds 



C cS = ^(V0) 2 + ^m 2 (A)0 2 (A) + ^A(A)0 4 (A) + 5C cS , 



(69) 
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where 5jC e g denotes all operators that are higher order in the and V0 and 
respect the symmetries. We have indicated that the parameters of the effective 
theory depend on the ultraviolet cutoff A. Sometimes we will omit this depen- 
dence for notational simplicity. One term that was omitted in (69) was /(A), 
which is the coefficient of the unit operator. This parameter can be interpreted 
as the contribution to the free energy from the scale T, or equivalently from 
the nonzero Matsubara frequencies. 

Including / in the effective Lagrangian, we can express the partition func- 
tion as a path integral over the three-dimensional field 0: 



: e ~f(A)V I V( p exp _ J d 3 x ^ 



(70) 



In the full theory, the partition function can also be written as a path integral: 



J V(f) exp - J dr J d 3 x C 



(71) 



Comparing (70) and (71), we obtain 
log Z = -f(A)V + logZ efr , 



(72) 



where Z e s is the partition function in the effective theory. In order to carry 
out matching, we treat the mass term in the effective theory as a perturbation. 
The Lagrangian is then divided according to 



(£eff) frec 1 
(£eff) int : 



2^ 



(73) 
(74) 



The mass parameter m 2 in (69) can be interpreted as the contribution to 
the Debye mass from the nonzero Matsubara modes or equivalently from the 
momentum scale T. The simplest way to determine the mass parameter is to 
calculate the Debye mass in the full theory and in the effective theory and 
demand that they be the same. In the original theory, the Debye mass is 
defined by the pole of the static propagator: 



p 2 + n(0,p) = 



V 



m 



D 5 



(75) 



where H(ou n , p) is the self-energy function. In the effective theory, the Debye 
mass is given by 



p 2 + m 2 + n cfr (p) = 



V 



m 2 D , 



(76) 



where n e ff (p) is the self-energy in the effective theory. The diagrams contribut- 
ing to the self-energy in the full theory are shown in Fig. 1 and the resulting 



22 



integral expression is given by 



Q Q 2 4 y TqrQIR* 

~6 g4 f QR Q 2 R 2 (P + Q + R) 2 ' (77) 

We have used the symbol ~ in (77) as a reminder that the relation only 
holds in strict perturbation theory. The last sum-integral depends on the ex- 
ternal momentum P, but it can be simplified by expanding it around p = 0. 
This expansion can be justified by noting that the leading-order self-energy 
is independent of the external momentum and that the leading-order Debye 
mass therefore is of order g 2 . To order g 4 , we can therefore set p = in the 
last two-loop sum-integrals. 

In the effective theory, the self-energy to two loops is given by the diagrams 
in Fig. 1 and is 

n cff ( P ) - / 1 - ^A 2 / 11 - ±A 2 / 11 1 + 5m 2 , (78) 
2 J q 2 4 J q 2 r 4 6 J q 2 r 2 p + q + r Y 

q qr qr 

where 5m 2 is a mass counterterm. The self-energy n e g(p) can similarly be 
expanded in powers of the external momentum p. To order A 2 , we can set 
p = 0. The external momentum p is the only scale in the loop integrals. 
When p = 0, these integrals vanish in dimensional regularization. Thus we 
can write 

IT cfr (0)^m 2 , (79) 
Matching the two expressions for the Debye mass, we obtain 



m 2 = 11(0,0) -5m 2 



-l 9 i QR Q 2 R 2 (Q + R) 2 - 6m2 - (80) 

Renormalization of the coupling constant g 2 is carried out by the substi- 
tution 



2(47r) 2 e ' 



^ = 1 + ^7- (81) 



After coupling constant renormalization, there is still a remaining pole in e. 
This divergence is cancelled by a mass counterterm which is 

Sm 2 = * , . 82 
24 4tt 2 e v ; 



23 



After mass renormalization, the mass parameter becomes 



^T 2 
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(83) 



The coupling constant g 2 satisfies a renormalization group equation: 



dg 2 Sg 4 



dn (4tt) ; 



(84) 



The logarithmic term — 3 log(/i/47rT) in (83) cancels the /x-dependence of the 
(yf 2 -term. The logarithmic term log(A/47rT) shows that the mass term de- 
pends explicitly on the ultraviolet cutoff A. This is necessary to cancel the 
A-dependence that arises in calculations using the effective theory. 

In the following, we need the coupling constant A in the effective theory 
only to leading order in the coupling constant g 2 in the full theory. At tree 
level we can simply read off A from the Lagrangian (1) of the full theory. 
Making the substitution $(t, x) — > y/T(f)(x.) and comparing Jq drC with the 
effective Lagrangian in (69), we obtain 



\ = 9 2 T. 



(85) 



We next consider the coefficient of the unit operator /. This parameter also 
has an expansion in powers of g 2 and can be interpreted as the contribution 
to the free energy logZ from the distance scale 1/T. It can be determined 
by calculating the free energy in the two theories and demanding that they 
be the same. In the full theory the Feynman diagrams contributing to log Z 
through three loops are shown in Fig. 5 (except for those with a cross): 
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The expression (86) is ultraviolet divergent. Renormalization is carried out 
by substituting the renormalization constant Z g 2 (81). The sum- integrals are 
given in appendix A and we obtain 
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(87) 
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In the effective theory, the diagrammatic expansion of the partition function 
is given by the diagrams in Fig. 5. Since the mass term is treated as a per- 
turbation, the resulting loop integrals are massless and therefore vanish in 
dimensional regularization. Eq (72) therefore reduces to 



V 

/(A) is then given by minus the left hand side of (87) divided by T. Using the 
renormalization group equation (84) it can be shown that /(A) is independent 
of the ultraviolet cutoff A up to corrections of order g 6 . 

4-1.2 Calculations in the effective theory 

Having determined the parameters in the effective Lagrangian, it can now 
be used to calculate physical quantities using perturbation theory. In order 
to take into account the physical effects of screening, we need to include the 
mass term in free part of the Lagrangian. The effective Lagrangian (69) is 
then decomposed as 

(^e ff ) frcc = i(V0) 2 + im 2 2 , (89) 

(A ff ) int = + 6C ■ (9°) 

Using this decomposition, we next consider the free energy. The diagrams 
contributing to the free energy through three loops are shown in Fig. 5 (except 
for those with a cross): 
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pqr v ' 

where 5m 2 is the mass counterterm given in (82). Using g 2 = XT together 
with the value for the mass counterterm and the results for the integrals listed 
in the appendix, we obtain 



logoff 1 3 A 2 y 
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V 12vr 8(4tt) 2 12(4tt) ; 



, A 9 , 

l0g 2^ + 8- l0g2 
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The free energy is the sum of (86) and (92). Substituting the expres- 
sion (83) for the mass parameter and expanding (92) in powers of g, we obtain 
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(93) 



Combining Eq. (86) and Eq. (93), we recover the weak-coupling expansion 
result (61) through order g 5 . This result is accurate up to corrections of order 
g 6 \ogg. Note that the explicit A-dependence of (92) cancels against the A- 
dependence of the mass parameter (83) to next-to-leading order in g. 

We have seen that the parameters /(A) and A(A) are independent of A, 
while m 2 (A) depends explicitly on the cutoff. The parameters of the effec- 
tive theory satisfy some evolution or renormalization group equations. These 
equations follow from the condition that physical quantities such as the Debye 
mass or the free energy are independent of the cutoff A. They can be derived 
by calculating the logarithmic ultraviolet divergences of the loop diagrams in 
the effective theory. For example, the coefficient of the unit operator satisfies 
the evolution equation 

dA ~ 12 (4tt)6 • (94j 

It can be derived as follows. At the four-loop level there are two diagrams 
contributing to the free energy that have logarithmic ultraviolet divergences. 

The poles of the first diagram are cancelled by the poles of the two-loop 
diagram with one mass counterterm insertion. The A dependence of the second 
diagram is cancelled by the A dependence of the parameter / and this yields 
the RG equation (94). Similarly, the mass parameter satisfies a renormalization 
group equation which is derived in the same manner and reads 

A lw=6(SF- (95) 

The solutions to (94) and (95) are trivial 

/(AWP.D-^log^, (96) 

m S(A) = m *(2,r) + 6Wlog4. (97) 

The solutions to the evolution equations can be used to sum up leading loga- 
rithms of g in physical quantities such as the free energy. This is accomplished 
by using the expressions (96) for the short-distance part of the free energy 
and (97) for the mass parameter of the effective theory: 
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(98) 



The first two terms in (98) is the contribution to the free energy from the 
scale T and include a term g 6 log g that arises from the solution (96) to the 
evolution equation. The other contribution to comes from expanding the Am 2 
term in (98) using (97). 

Fig. 6 clearly showed that the weak-coupling expansion is poorly behaved 
at large coupling. This fact motivated Blaizot, Iancu and Rebhan [130] to 
treat the soft sector differently. Instead of expanding the soft contribution to 
the pressure of QCD as a power series in g and truncating at the appropriate 
order, they calculated it in the loop expansion and kept the parameters of 
the effective Lagrangian as given by the dimensional reduction step. We apply 
the same strategy to the scalar field theory and show the resulting successive 
approximations in Fig. 8; at two loops we take the two- loop hard contribution 
and combine this with the two-loop soft contribution using the one-loop mass 
parameter. At three loops we take the three-loop hard contribution and com- 
bine this with the three-loop soft contribution using the two-loop mass. The 
bands are obtained by varying the renormalization scale fj, = A by a factor of 
2 around 2ttT. 

The three-loop pressure is then given by the sum of (86) and (92), where 
the parameters m 2 and A are given by (83) and (85). The result is shown in 
Fig. 8. It is interesting to note that the two-loop band is almost inside the 
three-loop band. The uncertainty in the three-loop approximation within this 
scalar toy model is quite large. Note that if the 0(g 6 log g) term is included 
as in Eq. (98) then the variation of the result with respect to the scale \i is 
less, namely, 0.93 < V/Videai < 0.95 at g — 4. However, in presenting the 
result shown in Fig. 8 we have not included this term because it formally 
contributes at four-loop order. Also note that within QCD the variation of 
the corresponding three-loop result is less owing to the independence of m E 
on the electric factorization scale, A E [130]. 

4.2 QED 

In this subsection, we will discuss dimensional reduction and the construc- 
tion of a three-dimensional effective field theory for QED. We call this effective 
field theory electrostatic QED (EQED) in accordance with the terminology in- 
troduced by Braaten and Nieto in the case of QCD [37]. The effective field 
theory is used to calculate the Debye mass and the free energy through or- 
der e 5 . In massless QED, there are two momentum scales that contribute to 
equilibrium quantities such as the free energy. In analogy with the scalar case, 
these are the scales T and eT, where e is the electromagnetic gauge coupling. 
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Fig. 8. The two- and three-loop approximations to the pressure where the soft part 
is not expanded out in powers of g. 

These scales are associated with the typical momentum of a particle in the 
plasma and the screening of static electric fields, respectively. 



4.2.1 Electrostatic QED 

The fields of the effective Lagrangian can, up to field redefinitions, be 
identified with the zero- frequency modes of the fields in QED. The timelike 
component of the gauge field, Aq, behaves as a real massive interacting field. 
The fact that A develops a thermal mass is a consequence of the breakdown 
of Lorentz invariance at finite temperature. Moreover, £ e ff is a gauge-invariant 
function of the spatial components A4. There is three-dimensional rotational 
symmetry and a discrete symmetry A — ► — A . The effective Lagrangian 
can then be written as 



C cS = X - (VA ) 2 + \m\A\ + \ E A* + X -Fl + 5C cS , (99) 



where <5£ c fr again denotes higher order operators that respect the symmetries. 
The operator Aq first contributes to the Debye mass at order e 5 and to the 
free energy at order e 6 . 

The coefficient of the unit operator /(A) was calculated in Ref. [35] through 
three loops: 



ll7r 2 T 4 j\ 50 20 

f (A) = < 1 a H 

J K ' 180 I 11 11 
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where a = e 2 /(4-7r) 2 . The mass parameter has been calculated by various 
authors through two loops [119], [35]: 



2 / a \ 167r2 ^2 

m E (A) = — - — T a 



i -I( iog 4^ + ^ +2iog2+ D 



a 



(101) 



The mass parameter coincides with the Debye mass through order a 2 . The 
reason is that there is no a 3//2 -term in weak-coupling expansion for the Debye 
mass. This is due to the facts that there are only fermionic propagators in the 
one-loop self-energy diagram in QED and that fermions need no resummation. 

The coupling constant is determined in the same manner as A in 4 - 
theory. The result is [106] 

e|(A) = e 2 T. (102) 

The coupling constant \e is calculated by evaluating the four-point function 
for timelike photons at zero external momenta. [106] Taking the different nor- 
malizations of the fields in full QED and EQED into account, one obtains 

W = -^. (103) 

Note that Xe requires no renormalization and vanishes at zero temperature. 
This reflects the gauge symmetry of QED. 

The renormalization group equation for the running gauge coupling is 

de 2 e 4 . 
"^=6^ < 104 > 

From Eq. (104), it is easy to check that the parameters /, m 2 E) e|, and \ E are 
in fact renormalization group invariant up to higher-order corrections. 

4.2.2 Calculations in EQED 

We next turn to calculations in the effective theory. In EQED, this is 
particularly simple since interactions first contribute to the free energy at 
order e 6 . Thus through e 5 , the contribution to the free energy is given by a 
simple one-loop calculation: 



log .Z, 
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i J \og(p 2 + m%) 

(105) 
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Expanding (105) in powers of e and truncating at order e 5 , we obtain 
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(106) 



The free energy through e 5 is the sum of Eqs. (100) and (106). This result 
was first obtained by Parwani using resummation. [34] In that calculation, the 
order e 5 contribution is obtained by expanding two- and three-loop diagrams 
as it was done in the scalar case in Sec. 3. Here, the result is obtained by 
a straightforward one-loop calculation and demonstrates that effective field 
theory significantly simplifies calculations of physical quantities. 

Similarly, the contribution to the Debye mass m 2 D from the scale gT is also 
given by a simple one-loop diagram in the effective theory. If the self energy 
in EQED is denoted by U c g(p, A), we have 



m 2 D = m 2 E + U eS (p,A) 



-m 2 E + 12\ E J 
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(107) 



The Debye mass to order e 5 was first obtained by Blaizot, Iancu, and Parwani 
using resummation techniques [119]. 

The resummation of logarithms of e is straightforward. The evolution 
equations for /(A) and m 2 E (A) are derived in the same manner as in the 
scalar case. The logarithmic divergences in the effective theory occur in the 
three-loop basketball diagram for / and the two-loop setting sun diagram for 
the mass parameter. These diagrams are both proportional to X 2 E . Thus the 
leading logarithmic term in the free energy is proportional to e 9 log e and in 
the Debye mass is proportional to e 8 loge. Thus the results for the free energy 
and Debye mass are correct through order O (a 3 log a). 



4.3 QCD 

In scalar field theory and QED, there are two momentum scales, T and 
gT, contributing to equilibrium properties. In nonabelian gauge theories, there 
are, in general, three momentum scales. They are 

• The scale T, which is the typical momentum of a particle in the plasma. 

• The scale gT, which is the scale associated with the screening of color- 
electric fields. 
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• The scale g 2 T , which is the scale associated with the screening of color- 
magnetic fields. 

The fact that there is an additional scale g 2 T in nonabelian gauge the- 
ories such as QCD, can be seen from their infrared properties. Perturbative 
QCD at high temperature is plagued by infrared divergences. The most se- 
vere divergences from a given infro ared divergent loop diagram comes from 
the region where all the Matsubara frequencies vanish and the loop momenta 
go to zero. Thus the infrared properties are those of pure-glue QCD in three 
spatial dimensions. This theory has a dimensionful coupling constant g 2 T and 
this is the only scale in the theory. It is a confining theory with a mass gap of 
order g 2 T that can only be calculated nonperturbatively. [39] 

This hierarchy of momentum scales suggests that one constructs a se- 
quence of two effective field theories that take care of the scales gT and g 2 T. 
This strategy was first proposed by Braaten and Nieto [37]. These effective 
fields theories are called electrostatic QCD (EQCD) and magneto static QCD 
(MQCD), respectively. The contribution to the free energy from the scale T 
is given by the coefficient of the unit operator /^(As) in EQCD. /^(As) has 
an expansion in powers of g 2 and starts to contribute to the free energy at 
order g°. The contribution from the scale gT is given by the coefficient of the 
unit operator /^(AgA^) m MQCD. t It has an expansion in powers of g and 
starts to contribute to the free energy at order g 3 . Finally, the contribution to 
the free energy from the scale g 2 T is given by the partition function in MQCD. 
It has an expansion in powers of g 2 , but first starts to contribute to the free 
energy at order g % '. 



4.3.1 Electrostatic QCD 

EQCD contains an electrostatic field Aq(x.) and a magnetostatic field 
Aj(x). The field Aq(x) behaves as a scalar field in the adjoint representation. 
The effective Lagrangian is 

£eqcd = \g%G% + ±(ViA*) 2 + \m\A%Al + \\ E {A%) 2 + 5C EQCD , (108) 

where G^ = diA a - — dj Af + g % f abc Af A a - is the magnetostatic field strength and 
DiAq = (<9j + gEt^A^AfyAQ is the covariant derivative. 5£eqcd consists of 
all higher-order gauge invariant operators than can be constructed out of the 
scalar field A and the gauge field Aj. In the matching procedure, the mass 
term in EQCD is again treated as a perturbation. 
The free energy for QCD can be written as 

F = Tf E (A E )- l ° gZ ° QCD (109) 

t The arbitrary factorization scales A# and Am separate the momentum scales T 
from gT, and gT from g 2 T, respectively. 
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where Zeqcd is the partition function for EQCD. The coefficient of the unit 
operator through three loops was calculated by Braaten and Nieto [37] 
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where a = g 2 /(An). A few comments are in order. After calculating the di- 
agrams that contribute to up to three-loop order and renormalizing the 
coupling constant by the substitution g 2 — ► Z 2 , where 

Z v~ l 3(4^ ' (1U) 

there is still a remaining pole in e. This pole is cancelled by the counterterm 
Sf E : 

Sf E = — • (112) 

l07T e 

This counterterm is analogous to the mass counterterm that was needed in 
the case of the scalar theory. The explicit A-dependence of f(A E ) is necessary 
to cancel logarithms of A/itle arising in calculations of vacuum diagrams in 
EQCD. The coefficient of the unit operator therefore satisfies the evolution 
equation 

d dACA 2 / 110 \ 

AE dT E fE = -— amE - (113) 

Note also that we need the order e of the mass parameter m 2 E since it multiplies 
the pole in e in (112). The mass parameter vtie is calculated by matching the 
Debye mass in QCD and EQCD using strict perturbation theory. * The result 
through next-to-leading order has been calculated by Braaten and Nieto [37] 
through order e. The order e° term is 



* Note that the Debye mass cannot be defined perturbatively beyond leading order 
in g due to infrared divergences associated with the magnetostatic scale g 2 T [113] 
and therefore requires a nonperturbative definition [105]. See also Sec. 7. 
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The order-e term is 
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The coupling constants qe and A# are calculated in the same manner as in 
QED. The result is 
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The gauge coupling g 2 satisfies the RG-equation 



(116) 
(117) 
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118) 



Using (118), it can be shown that the mass parameter m 2 E as well as the 
coupling constants g\ and \ E are independent of the factorization scale A E - 
Having determined the parameters in EQCD, one can calculate the con- 
tribution to various physical quantities from the scale gT. In order to do so, 
we again include the effects of the mass term for Aq into the free part of the 
Lagrangian of EQCD. The contribution to the free energy from the scale gT is 
given by / which is equal to the logarithm of the partition function of EQCD: 



M 



T log Z EQCD 
V 



(119) 



The calculation of ^eqcd through three loops was done by Braaten and Ni- 
eto. [37] The result is 
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(120) 
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By expanding (120) in powers of g E , substituting (114) and (116), we obtain 
the contribution to the free energy from the scale gT through g 5 : 
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(121) 



where a s = (C A + Tp)a s /37r. The free energy through g is the sum of (110) 
and (121) and is given by (65) for C A — N c — 3. Note that the dependence on 
factorization scale A E cancels between f E and f M up to corrections of order 

9 6 - 



4.3.2 Magnetostatic QCD 

The effective Lagrangian of magnetostatic QCD is 



-MQCD 



MQCD , 



(122) 



where = diA a - — djA": + g M f ahc A^A^ is the magnetostatic field strength 
and gM is the gauge coupling. The term <5£mqcd consists of all possible gauge- 
invariant terms that can be constructed out of the spatial gauge field Aj. In 
Ref. [37] it was shown that the operators in 5C first contribute to the free 
energy at order g 12 . In the remainder of this section, we therefore omit 5C and 
only consider the minimal gauge theory with action / d 3 x \G^G^. 

For the present discussion, the only parameter we need to determine is 
the gauge coupling. To leading order in g E , one finds [37] 



gM — 9e 



(123) 



To this order in g, the parameter gM is independent of the factorization scale 
Am- The free energy of QCD can be written as 



F = Tf E (A E )+Tf M (A E ,A M ) 



T 'log Z MQCD 

V 



(124) 



where /ju"(A £ , A M ) is the unit operator of MQCD and logi^g^ is the par- 
tition function in magnetostatic QCD. 
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Perturbation theory using the Lagrangian (122) is plagued with infrared 
divergences. In the case of the free energy, these infrared divergences first 
enter at the four-loop level. Thus perturbation theory can not be used to 
calculate the contribution to the free energy from the scale g 2 T. However, for 
dimensional reasons, the nonperturbative contribution has to be of the form: 

\ogZ MQCD ( A M \ 6 

77 = U + Mog— \g b M , (125) 

v \ 9 m J 

where a and b are pure numbers. The coefficient a can be determined by 
using nonperturbative methods such as lattice simulations. The coefficient 
b can be calculated by extracting the logarithmic ultraviolet divergences in 
the four-loop diagrams contributing to the free energy in MQCD. This was 
recently done in a calculation by Kajantie, Laine, Rummukainen and Schroder 
[107,108]: 

h J r* 9m /43 157vr 2 \ 

h = dACA Wr\T2-^m) ■ (126) 

The free energy to order g 5 has been calculated both by resummation and 
effective field theory methods. Using effective field theory, one can also deter- 
mine the coefficient of g 6 logg. There are two contributions to this coefficient. 
The first contribution comes from the scale gT and can be determined by 
perturbative calculations of the logarithm of the partition function in EQCD. 
It involves the evaluation the of four-loop vacuum diagrams and has recently 
been determined in a very impressive calculation by Kajantie, Laine, Rum- 
mukainen and Schroder: [107,108] 

f (4)/ A * \ j /ht3 9 m (43 491tt 2 \ A E 

f M (A*, Am) - d A C AJ ^- A - —j log — , (127) 

where the superscript indicates the number of loops. The second contribution 
comes from the scale g 2 T and is given by the coefficient b in Eq. (126). 



5 Screened Perturbation Theory 

In this section, we discuss one possibility to reorganize perturbation theory, 
which is screened perturbation theory (SPT) which was introduced by Karsch, 
Patkos and Petreczky [86] . It can be made more systematic by using a frame- 
work called "optimized perturbation theory" that Chiku and Hatsuda [87] 
have applied to a spontaneously broken scalar field theory. In this approach, 
one introduces a single variational parameter m, which has a simple interpre- 
tation as a thermal mass. The advantage of SPT is that it is relatively easy 
to apply and that higher-order corrections can be calculated so that one can 
study convergence properties. 
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Fig. 9. Diagrams which contribute up to three- loop order in screened perturbation 
theory. A boldfaced x indicates an insertion of m 2 . 

In SPT a mass term is added and subtracted to the scalar Lagrangian 
with the added bit kept as part of the free Lagrangian and the subtracted 
bit associated with the interactions. This can be accomplished by introducing 
a parameter S in the spirit of the linear delta expansion [85]. The resulting 
Euclidean Lagrangian is 



where we have also introduced a factor of 5 in the <ft A term. In the limit 5 — > 1 
the SPT Lagrangian reduces to the standard scalar Lagrangian. 

The Lagrangian (128) is then expanded into free and interacting parts 
with all interaction terms proportional to 5 



The loop expansion then simply becomes an expansion in S. Note also 
that the counterterms that are necessary to renormalize the theory also have 
expansions in powers of 5. After expansion of all of the contributions to the 
appropriate order in 5 we set 5 = 1 . If we were able to calculate the result to all 
orders in 5 the result would independent of the mass parameter m; however, 
at any finite order in the expansion there will be a residual dependence on 
m. We will discuss three possible ways to fix the parameter m including a 
variational prescription. The diagrams which contribute up to the free energy 
up to three-loop order in SPT are shown in Fig. 9. 

5.1 One-loop contribution 
The one-loop free energy is 



Using the results from Ref. [88] for the sum-integrals and counterterm 
A £o 5 we obtain the final result for the one- loop free energy 



(128) 




(129) 
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(47t) 2 ^ = -\(2L + 3)m 4 - ± J T 4 , (132) 

where J„ are functions of m/T defined in (A. 10) evaluated at e = and 
L = log(jU 2 /m 2 ). 



5.2 Two-loop contribution 

The contribution to the free energy of order 5 is 



f)7- 

Tx = Ti, + Fib + A^o + Aim 2 , 



(133) 



where Ai£o and Aim 2 are the order-5 vacuum and mass counterterms, respec- 
tively. Using the results from Ref. [88] for the sum-integrals and counterterms, 
the two-loop contribution to the free energy becomes 



(47r) 2 J- 1 = -[(L + l)m 2 - JxT 2 
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(134) 



5.3 Three-Loop contribution 



The contribution to the free energy of order S 2 is 
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where we have included all of the order-5 2 counterterms. Using the results 
from Ref. [88] for the integrals and counterterms the three-loop contribution 
to the free energy becomes 
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+ (3(3L + 4)J 2 + 3J 2 J 2 + 6ir 2 + 4K3)T 4 ] , (136) 



where C = 39.429, C\ = —9.8424, and K 2 and K 3 are functions of m/T given 
in Ref. [120]. 
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5.4 Pressure to three loops 



The contributions to the free energy of zeroth, first, and second order in 5 
are given in (132), (134), and (136), respectively. Adding them we obtain the 
successive approximations to the free energy in screened perturbation theory. 
Using the fact that V = —T we obtain the one-loop approximation to the 
pressure: 



(4tt) 2 P = o [ 4 J oT 4 + (2L + 3) m 4 

8 L 

The two-loop approximation is obtained by including (134): 
1 
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The three-loop approximation is obtained by including (136) to obtain: 
1 
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+ (5L 3 + 17L 2 + f L - 23 - ffvr 2 - <(1) + C ) m 

(139) 

The only remaining task is to determine the mass parameter m which we will 
discuss next. 



5.5 Mass prescriptions 

The mass parameter m in screened perturbation theory is completely ar- 
bitrary. In order to complete a calculation using SPT, we need a prescription 
for the mass parameter m = m*(T). The prescription of Karsch, Patkos, and 
Petreczky for (T) is the solution to the one-loop gap equation: 



J 1 (tti,/T)T 2 



21 g^- + l)m: 



(140) 



where the function Ji is defined in (A. 10). Their choice for the scale was 
/i* = T. In the weak-coupling limit, the solution to (140) is m* = g(/i*)T '/ \^2A. 
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There are many possibilities for generalizing (140) to higher orders in g. We 
will consider three different possibilities in the following. 



5.5.1 Debye mass 

One class of possibilities is to identify m* with some physical mass in the 
system. The simplest choice is the Debye mass m D defined by the location of 
the pole in the static propagator: 



The Debye mass is a well defined quantity in scalar field theory and abelian 
gauge theories at any order in perturbation theory. However, in nonabelian 
gauge theories, it is plagued by infrared divergences beyond leading order. [113] 

5.5.2 Tadpole mass 

The tadpole mass is another generalization of Eq. (140) to higher loops. 
The tadpole mass is defined by 



The tadpole mass is well defined at all orders in scalar field theory, but the 
generalization to gauge theories is problematic. The natural replacement of 
((p 2 ) would be (AftAfj), which is a gauge- variant quantity. 

5.5.3 Variational mass 

There is another class of prescriptions that is variational in spirit. The 
results of SPT would be independent of m if they were calculated to all or- 
ders. This suggests choosing m to minimize the dependence of some physical 
quantity on m. The variational mass is defined by minimizing the free energy; 



The variational mass has the benefit that it is well defined at all orders in 
perturbation theory and can easily be generalized to gauge theories. 

5.5.4 Comparison 

At one loop, the three different prescriptions give the same gap equation, 
Eq. (140). Moreover, it turns out that the two-loop tadpole mass coincides with 
the one-loop tadpole mass. However, at two-loops the screening and variational 
masses are ill-behaved [88] . The screening mass solution ceases to exist beyond 
g ~ 2.6 and the variational gap equation has no solution for finite g. In Fig. 10 
we show the various solutions m* (T) to the SPT one- loop gap equation (shaded 
band) and two-loop gap equations (lines) as functions of g(2nT). From this we 



p 2 + m 2 + S(0,p) = , 



2 2 
p = —m D . 



(141) 




(142) 
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Fig. 10. Solutions m*(T) to the SPT one-loop gap equation (shaded bands) and 
two- loop gap equations (lines) as functions of g(2irT). 

conclude that the best prescription is to use the tadpole mass which results 
in the same gap equation at one- and two-loops given by Eq. (140). 

5. 6 Pressure 

In Fig. 11, we show the one-, two-, and three- loop SPT-improved ap- 
proximations to the pressure using the tadpole gap equation. The bands are 
obtained by varying fi by a factor of two around the central values \i = 2nT. 
The one-loop bands in Fig. 11 lie below the other bands; however, the two- 
and three-loop bands both lie within the g 5 band of the weak-coupling expan- 
sion in Fig. 6. The one-, two-, and three- loop approximations to the pressure 
are perturbatively correct up to order g 1 , g 3 , and g 5 , respectively; however, 
we see a dramatic improvement in the apparent convergence compared to the 
weak-coupling expansion. Improved convergence properties were also found 
for the entropy and the Debye mass [88]. 

This result demonstrates the effectiveness of SPT in providing stable and 
apparently converging predictions for the thermodynamic functions of a mass- 
less scalar field theory. An essential ingredient of this approach is using the 
solution to a gap equation as the prescription for the mass parameter m. In the 
following section we will see how this idea can be extended to gauge theories. 

6 Hard-thermal-loop Perturbation Theory 

In this section, we discuss the generalization of SPT to gauge theories. 
This framework is called hard-thermal- loop perturbation theory (HTLPT). 




40 



7iT < n < 471T 




0.9 



1 2 

g(2;iT) 



Fig. 11. One-, two-, and three-loop SPT-improved pressure as a function of g(2irT). 
Bands correspond to a variation of the renormalization scale irT < fj, < AttT. 

The HTLPT framework is attractive because it allows for a systematic ana- 
lytic reorganization of perturbative series based on the HTL effective action. 
Additionally, it is manifestly gauge invariant and is applicable to calculating 
both static and dynamical quantities. For simplicity, here we restrict the dis- 
cussion to pure-glue QCD. Calculations including fermions can be found in 
Refs. [91,93]. 

The Minkowski-space Lagrangian density that generates the perturbative 
expansion for pure-glue QCD can be expressed in the form 

£ Q cd = "^Tr {G^GT) + £gf + £ ghost + A£ QCD , (144) 

where G^ = d fl A u — d u A fl — ig[A f j i , A u ] is the gluon field strength and A^ is the 
gluon field expressed as a matrix in the SU (N c ) algebra. The ghost term £ g h os t 
depends on the choice of the gauge-fixing term £ gt . In HTLPT the Lagrangian 
density is written as 

C = (£qcd + £htl) rx + AAttl • (145) 

The HTL improvement term appearing above is (1 — 5) times the isotropic 
HTL effective action which generates all HTL n-point functions[109] 

£htl = - 5)m 2 D TT (g^ (jpj^ G%J , (146) 

where is the covariant derivative in the adjoint representation, = (1, y) 
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is a light-like four- vector, and (. . .) y represents the average over the directions 
of y. The free Lagrangian in general covariant gauge is obtained by setting 
5 = in (145): 



C ircc = -Tr (d^A u d^A u - d,A u d»A») - -Tr [(d^f 



— -m 2 D Tr 



(d^A a — daAfj) 



y a y (5 
(ydf 



(d»Ap - dpA^ 



(147) 



The resulting propagator is the HTL gluon propagator and the remaining 
terms in (145) are treated as perturbations. The propagator can be decom- 
posed into transverse and longitudinal pieces which in Minkowski space are 
given by 



A T (p) 
Al(p) 



1 



p 2 - n T (p) ' 
i 

-nip 2 + n £ (p) ' 



(148) 
(149) 



where TIt and 11^ are the transverse and longitudinal self-energies, respec- 
tively, and the four- vector is 



U U n 'P 1 



(150) 



The four- vector n M specifies the thermal rest frame (canonically n — (1,0)). 
In terms of II T and Yi L the full self-energy tensor is given by 



1 



n^(p) = -U T (p)T p ^ - -U L (p)L^ , 



(151) 



where the tensors T p and L p are 



' V ni 



n v 



(152) 
(153) 



The HTL gluon self-energy tensor for a gluon of momentum p is given in 
Eq. (22) with m 2 D given by Eq. (31) and in the limit that d — > 3 the resulting 
expressions for the transverse and longitudinal gluon self-energies reduce to the 
expressions given in Eqs. (29) and (30) with m 2 D given by Eq. (32). Note that 
there are also HTL vertex corrections which are given by similar but somewhat 
more complicated expressions which can be found in Refs. [28,92,93]. 
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Fig. 12. Longitudinal and transverse dispersion relations. 



As mentioned above, HTLPT is a systematic framework for performing 
calculations in thermal gauge theories which is gauge invariant by construc- 
tion. It systematically includes several physical effects of the plasma such as 
the propagation of massive particles, screening of interactions, and Landau 
damping. We briefly discuss these next. 

Massive quasi-particles: 

The HTL self-energies are included in the zeroth order propagators which 
results in the resummed HTL propagator for gluons. The dispersion relations 
for the transverse and longitudinal gluonic degrees of freedom are determined 
by locating the zeros of the inverse propagator which gives the following two 
equations 



The dispersion relations for transverse and longitudinal gluons are shown in 
Fig. 12. As can been seen from this Figure both modes approach a constant 
in the limit of small momentum and approach the light-cone in the limit of 
large momentum. 

Screening: 

HTLPT also includes screening of interactions which can be seen by ex- 
amining the static limit of the HTL propagators. For instance, the inclusion 
of the longitudinal self-energy changes the Coulomb potential of two static 
charges in the plasma to a Yukawa potential: 



— k 2 — Ut(ut, k) — 0, 
k 2 + U L {u L ,k) = . 



(154) 
(155) 



lim Al(u, k) 



1 



(156) 



k 2 + m 



2 ' 
D 
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This result shows that chromoelectric fields are screened on a scale r ~ vtv^ . 
Likewise, the screening of long wavelength chromomagnetic fields is deter- 
mined by the transverse propagator for small frequencies 



This shows that there is no screening of static magnetic fields meaning that 
the magnetic mass problem of nonabelian gauge theories at high temperature 
is not solved by HTL resummation. However, HTL resummation does give 
dynamical screening at a scale r ~ (m 2 D uj) 3 . Note that the divergences asso- 
ciated with the absence of static magnetic screening do not pose a problem 
until four-loop order; however, at four loops the lack of static magnetic screen- 
ing gives rise to infrared divergences that cause perturbation theory to break 
down. 

Landau damping: 

The transverse and longitudinal HTL self-energies also contain the physics 
of Landau damping. Landau damping represents a transfer of energy from the 
soft modes to the hard modes for spacelike momentum. This can be seen 
from the analytic structure of the self-energies given by Eqs. (29) and (30). 
Because of the logarithms appearing in these functions there is an imaginary 
contribution to the self-energies for —k<u<k which gives the rate of energy 
transfer from the soft to hard modes. Note that ignoring this contribution leads 
to gauge variant and unrenormalizable results. 

6.1 Calculation of the free energy using HTLPT 

We will now present results of a two-loop calculation of the free energy of 
pure-glue QCD within HTLPT. The relevant one- and two-loop diagrams are 
shown in Fig. 13 with shaded circles indicating dressing of propagators/ vertices 
with HTL resummed propagators/vertices. 



Fig. 13. Diagrams contributing through NLO in HTLPT. Shaded circles indicate 
dressed HTL propagators and vertices. 

6.1.1 One-loop contribution 

The thermodynamic potential at leading order in HTLPT is 




(157) 




{N 2 C - \)T g + A S , 



(158) 
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where Ao£o is a vacuum counterterm and T g is shown in Fig. 13. Using the 
results from Ref. [92] for the sum-integrals and counterterms the one-loop free 
energy becomes 

^lo , 15 „ 2 , on - 3 
- 1 m n + 30m, 



v - ■ l D ^ UU '"D 

ideal ^ 

+ f (; + 21og§-7 + 2 7 + ^m 4 D , (159) 

where m = mi)/2iiT and jx = /z/27rX '. 

6.1.2 Two-loop contribution 

The diagrams that contribute to the two-loop free energy are shown in 
Fig. 13. The thermodynamic potential at next-to-leading order in HTLPT 
can be written as 

^nlo = ^lo + (N 2 C - 1) [F 3g + T Ag + T gh + F ct ] + A^ 

+A 1 m^^- r fi LO , (160) 

where A]£o and &.\m 2 D are the terms of order S in the vacuum energy density 
and mass counterterms, respectively. Again using the results from Ref. [92] 
for the sum-integrals and counterterms the two-loop contribution to the free 
energy can be computed. Combining this with the one-loop contribution gives 
the total next-to-leading order result for the thermodynamic potential 



^NLO 1 . _ . 3 45 / fl 7 7T 2 

— = 1 — 15mn log h 7 H 

^idcai ° 4 I S 2 2 3 



N c a s 
3tt 



15 Jr . 165 / a 36, . n n \ „ , 
- — + 45m D - — I log | - — logm D -2.001 I m 2 D 




log 2 + 22+7)^ 



(161) 



In Fig. 14, we have plotted the LO and NLO HTLPT predictions for the 
pressure of pure-glue QCD as a function of T/T c , where T c is the critical 
temperature for the deconfinement phase transition. We have also included 
the 4-dimensional lattice gauge theory results of Boyd et al. [51] and the 3- 
dimensional lattice gauge theory results of Kajantie et al. [95] in this Figure. 

To translate a s (27rT) into a value of T/T c , we use the two- loop running 
formula for pure-glue QCD with = 1.04 T c . Thus a s {2nT) = 0.06 and 
0.2 translate into T/T c ~ 400 and T/T c ~ 1, respectively. The LO and NLO 
HTLPT results are shown in Fig. 14 as a light-shaded band outlined by a 
dashed line and a dark-shaded band outlined by a solid line, respectively. The 
LO and NLO bands overlap all the way down to T = T c , and the bands 
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Fig. 14. The LO and NLO results for the pressure in HTLPT compared with 4-d 
lattice results (diamonds) and 3-d lattice results (dotted lines) for various values of 
an unknown coefficient in the 3-d effective Lagrangian. The LO HTLPT result is 
shown as a light-shaded band outlined by a dashed line. The NLO HTLPT result is 
shown as a dark-shaded band outlined by a solid line. The shaded bands correspond 
to variations of the renormalization scale /j, by a factor of two around [i = 2itT. 

are very narrow compared to the corresponding bands for the weak-coupling 
predictions in Fig. 7. Thus the convergence of HTLPT seems to be dramati- 
cally improved over naive perturbation theory and the final result is extremely 
insensitive to the scale fi. 

Despite the apparently improved convergence, the predictions of HTLPT 
seem to lie significantly above the available 4d lattice results available for 
T < 5T C . This could be due to at least three possibilities: (a) the expansion 
in the variational mass performed is poorly convergent; (b) by including the 
HTL resummation also for the hard modes HTLPT treats the hard modes 
improperly suggesting that a method which explicitly treats these separately 
would be better; or (c) HTLPT is missing crucial physics related to confine- 
ment near T ~ T c . Out of these three possibilities the last two are the most 
plausible explanations given the success of such mass expansions applied to 
scalar theories [89]. We will discuss an alternative method of reorganization 
which addresses the second possibility in Sees. 8 and 9. The third possibility 
could, in principle, be taken into account by using an effective theory which is 
written in terms of Polyakov loops [121]. For a particularly nice introduction 
see Ref. [122]. Unfortunately, work on this is still unable to produce a method 
for calculating the coefficients in the Polyakov loop effective Lagrangian near 
T c where it would be needed. Note though that there have been recent ad- 
vances in the calculation of the one-loop Polyakov effective action [123,124] 
and attempts to extract renormalized Polyakov loops from the lattice [125]. 
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Also there has been an attempt to introduce models which include a non-trivial 
Polyakov loop expectation value in a quasiparticle-based model [126]. 

The reader may also notice that the HTLPT results also seem to disagree 
with the 3d lattice calculations shown in Fig. 14. In this context it is important 
to note that the lattice results from Ref. [95] are now known to be incorrect due 
to the assumption that one of the coefficients in the 3d effective Lagrangian, 
namely /3ei, was zero. It is now known to be non-zero and this is expected to 
modify the result significantly at low temperatures [108]. 



7 ^-derivable Approach 

In the previous sections, we have discussed SPT and HTLPT, which are 
approaches that involve the variational determination of a mass parameter 
m 2 . In this section, we will discuss the ^-derivable approach, which is another 
approach that is variational of nature. In this approach, one uses the exact 
propagator as a variational function. Its formulation was first constructed by 
Luttinger and Ward [67] and by Baym [68]. Later it was generalized to rela- 
tivistic field theories by Cornwall, Jackiw and Tomboulis [69]. The approach 
is based on the fact that the thermodynamic potential can be expressed in 
terms of the 2 particle-irreducible (2PI) effective action which has a diagram- 
matic expansion involving the 2PI skeleton graphs. Although here we focus 
on equilbrium physics we also note that the 2PI formalism and its general- 
izations are also very useful when studying non-equilibrium real-time physics 
[127,128,129]. 

Applying the ^-derivable approach to quantum field theories, one is fac- 
ing two nontrivial issues. The first issue is the question of renormalization. 
The three-loop calculation by Braaten and Petitgirard [80] for scalar field the- 
ory shows that there are ultraviolet divergences at order g 6 that cannot be 
eliminated by renormalization. This explicit calculation seems to be in con- 
tradiction with the papers by van Hees and Knoll [76] and by Blaizot, Iancu 
and Reinosa [77] , which show that ^-derivable approximations can be system- 
atically renormalized. 

The second issue is that of gauge dependence. While the exact 2PI effective 
action is gauge independent at the stationary point, this property may be lost 
in approximations. This problem has recently been studied by Arrizabalaga 
and Smit [78], where they showed that the n-loop ^-derivable approximation 
<£>„, defined as the truncation of the action functional after n loops, has a 
controlled gauge dependence. The gauge dependence of $ n shows up at order 
g 2n , where g is the gauge coupling. Moreover, if the nth order solution to 
the gap equation is used to evaluate the complete effective action, the gauge 
dependence first shows up at order g 4n . 

The ^-derivable approach has several attractive features. First, it respects 
the global symmetries of the theory. Thus it is consistent with the conservation 
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laws that follow from the Noether's theorem. Second, when evaluated at the 
stationary point, one is guaranteed thermodynamic consistency. [68] Moreover, 
it turns out that the two-loop ^-derivable approximation has an additional 
property. The entropy reduces to the one-loop expression at the variational 
point. This property was first shown for QED by Vanderheyden and Baym 
[110], and later generalized to QCD by Blaizot, Iancu and Rebhan [72]. 

Since the ^-derivable approach involves the propagator as a variational 
function and not simply a mass parameter, it is also difficult to solve. Braaten 
and Petitgirard [80] have therefore developed a strategy to systematically solve 
the n-loop ^-derivable approximation to the thermodynamics of massless 4> 4 
field theory. The method involves expanding sum- integrals in powers of m/T 
and g 2 , where m is a variational mass parameter. Blaizot, Iancu and Reb- 
han [72] have solved approximately the two-loop ^-derivable approximation 
in QCD by exploiting the fact that the hard-thermal-loop self-energies are the 
solutions to the gap equations to leading order in g for soft external momenta. 
We discuss these approximations next. 

7. 1 Scalar field theory 

We start our discussion by again considering a massless thermal scalar 
field theory with a Lagrangian given by (1). The thermodynamic potential is 



where D is the exact propagator, II is the exact self-energy, and 3>[D] is the 
sum of all two-particle irreducible diagrams. The two-particle irreducible di- 
agrams are shown diagrammatically in Fig. 15 up to three-loop order. The 




i 




(162) 




Fig. 15. ^-derivable two- and three-loop skeleton graphs. 



exact propagator D(P) satisfies the variational equation 



5Q[D] 
5D 



. 



(163) 



Using (162) the variational equation (163) can be written as 



n(P) = 2 



6$[D] 
5D{P) ' 



(164) 



In the case of a thermal scalar field, the exact propagator is 
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(P) = 

1 ; P 2 + n(P) 



(165) 



Substituting (165) into (162), the thermodynamic potential Q[D] becomes 

m = \f p \og (P 2 + H(P)) - \i p n{P) p2 ^ {p) + *[D] . (166) 

If we truncate $ at the nth order in the loop expansion, one refers to the corre- 
sponding variational approximation fl n to the exact thermodynamic potential 
Q as the n-loop ^-derivable approximation. The n-loop ^-derivable approx- 
imations, where n = 2,3,... define a systematically improvable sequence of 
approximations to the exact thermodynamic potential. 



7.1.1 Two-loop ^-derivable approximation 

In the two-loop ^-derivable approximation there is only a single diagram 
contributing to $[D] which is the first diagram on the right-hand side in 
Fig. 15. The thermodynamic potential f2 2 is 



fi 2 [^] = ^ p iog(p 2 + n(P))-^ 



IT(P) 



P p 2 + n(P) 

v 2 fc^W)) 2 - (167) 

It follows from Eq. (164) that the self-energy is independent of the momentum 
and we shall denote the exact self-energy IT in this approximation by m 2 . The 
one-loop gap equation becomes 

m 2 = \g 2 i D2 1 2 ■ (168) 

2 TpP 2 + m 2 

One can use the one-loop gap equation to simplify the expression for the 
thermodynamic potential: 

This eliminates all reference to terms which are explicitly of order a. 

The two sum-integrals are given in (A. 7) and (A. 8). They are separately 
ultraviolet divergent, but the sum is finite. We can take the limit e — > and 
obtain 



1 



8(4tt) 



4 J T 4 + 2 J im 2 T 2 + m 4 . (170) 



This result is identical to the free energy in SPT at two loops. The gap equa- 
tion (168) is ultraviolet divergent and requires renormalization. We write the 



renormalization constant Z g 2 as 
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Z 9 2g 2 = g 2 + A ig 2 + A 2 g 2 + ... , (171) 

where the subscript indicates the loop order. Making the substitution g 2 — > 
Z g 2g 2 , the gap equation (168) can be written as 




" 2TP 2 + m 2 ' 
Using Eq. (A. 9) in the appendix and 



(172) 



the gap equation becomes finite: 



2 1 

m = -a 
2 



JiT 2 -(L + l)m 2 . (174) 



Here L = log(/i 2 /m 2 ). 

It is interesting to note that both the renormalized thermodynamic po- 
tential and renormalized gap equation are identical to the expressions in 
SPT, Eq. (134) and (140). However, the coupling constant in the two-loop 
^-derivable approximation runs differently with the renormalization scale \x: 

^ = Wf (175) 

Hence, the /3-function differs from the perturbative one (62) by a factor of 
three. 

The two-loop ^-derivable approximation to the entropy can be found by 
differentiating the thermodynamic potential (167) with respect to the tem- 
perature and using the gap equation to eliminate the ultraviolet divergences. 
This gives 



TS 



1 



(4tt) ; 



2J T 4 + Jim 2 T 2 . (176) 



This expression is equal to the entropy of an ideal gas of massive bosons. 

The fact that the two-loop expression for the entropy reduces to the one- 
loop expression does not depend on the fact that the self-energy is a constant 
mass term. This simplification at the two-loop level motivated Blaizot, Iancu 
and Rebhan [72] to use the expression for the entropy as the starting point 
for various approximations to the thermodynamic functions. In the case of a 
scalar field theory, we have 

s = -§rf P log ( p2 + ™ 2 ) > ( 177 ) 
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where m 2 again is a thermal mass. If we use the weak-coupling result for 
m 2 , Eq. (177) reproduces the weak-coupling result for the entropy through 
order g 2 . However, it underestimates the order- g 3 correction to the entropy 
by a factor of four [72] . The reason is that the weak-coupling result for m 2 is 
the solution to the gap equation only to order g 2 . Blaizot, Iancu and Rebhan 
corrected for that by using the order-g 3 solution to the gap equation: 



2tt 2 



m 



-aT 2 



1 - V&a 1 ' 2 



(178) 



The approximation using Eq. (178) for the mass obviously breaks down 
for sufficiently high values of g, since the mass turns negative. The authors of 
Ref. [72] then introduced the mass 



2 2tt 2 2 
m = ——al 



VGa 1 / 2 



(179) 



in their approximation to entropy. This approximation gives the same para- 
metric improvement as (178). Finally, they introduced what they called the 
next-to-leading order HTL approximation by truncating the gap equation af- 
ter terms that are third order in g and m/T and solving the resulting quadratic 
equation with respect to m. This amounts to using the mass parameter 



m = 



V6 



na 1/2 Jl + 3a/327r 2 



1 



-a 



T 



;iso) 



The entropy obtained by using the mass prescriptions (179) or (180) is non- 
perturbative in the sense that there are contributions from all orders in g. 



7.1.2 Three-loop ^-derivable approximation 

The three-loop ^-derivable approximation to the free energy is 




vp 2 + it(p) j 

1 4 r 1 1 1 1 

~48 9 TpqrP 2 + n(P) Q 2 + IT(Q) R 2 + U(R) S 2 + U(S) ' (181) 

where S — —(P + Q + R). The corresponding gap equation that follows from 
the variation of (181) with respect to n(P) is 



1 



1 



j qQ 2 + IT(Q) 6 y TqrQ 2 + U(Q) R 2 + U(R) S 2 + U(S) ' 

(182) 
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It follows from the gap equation (182) that the self-energy II(P) is a nontrivial 
function of the external momentum P. 

In the previous subsection, we saw that we could use the gap equation to 
eliminate the ultraviolet divergences in two-loop approximation to the ther- 
modynamic potential. We will follow that strategy again and use (182) to 
eliminate some of the divergences in (181). We then obtain 



^ = ^ p iog(p 2 + n(P))-^ 



n(p) 



pP 2 + n(P) 

+ 48^pq;?p 2 + n(P) q 2 + n(Q) r 2 + u(R) s 2 + u(S) ' (183) 

In section 5, we saw that a reasonable approximation in SPT is to expand the 
sum-integrals in powers of m/T treating m to be of order gT. A similar strat- 
egy was developed and applied to the three-loop ^-derivable approximation by 
Braaten and Petitgirard [80]. The idea is to introduce a mass variable which 
in the weak-coupling limit is of order gT, and then expand the sum-integrals 
as double expansions in g 2 and m/T . Braaten and Petitgirard chose to use the 
Debye mass as a mass variable. In the three-loop ^-derivable approximation, 
the Debye mass is the solution to the gap equation 

m2 =h%Q^IHQ)-l^'' (184) 
where 

2bebye = ^__^ 9! , TT77T\ 7n , -rr/m , tt/ca ■ ( 185 ) 



Q rQ 2 + n(Q) r 2 + n(p) s 2 + u(S) 



p=im 



In the three-loop ^-derivable approximation, the variational equation (182) 
can be rewritten as 

n(P) = m 2 - ^g A [J sun - J Dcbyc ] . (186) 
where 

Xs ^ P) r QR Q 2 + n(Q) r 2 + n(p) s 2 + u(s) ' (187) 

Since we are assuming that the scales m ~ gT and 27rT can be separated, 
we are allowed to expand II(P) in powers of g and m/T. The variational gap 
equation will be solved in the two momentum regions separately. For hard 
momenta, the self-energy is expanded as 

n(p) = m 2 + g " [n 4 , (P) + n 4 ,i(P) + ...] + ..., (iss) 
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where the functions II nj fc(P) are of T 2 (m/T) k when the external momentum P 
is hard. For soft momentum P = (0, p), the self-energy is written as 11(0, p) = 
m 2 + cr(p), where 

*(p) =9 4 K- 2 (p) + <J±M + ...]• (189) 

Here, the functions u n ^{p) are of order m 2 (m/T) k when p is soft. By in- 
serting the expansions for II(P) and a(p) into the gap equation and expanding 
systematically in g and m/T, we can find expressions for H n ^{p) and er n ,fc(p) 
by matching coefficients of order g n on both sides and solving it recursively. 
In this manner, we find 

\r i i 

IMP) = -g} qrQ 2 R 2 {P + Q + R) 2 + 6 T2/DC ^ ' ^ 
<7 4> - 2 (p) = [/sun " W] , (192) 

where J sun and ioebye are the three-dimensional counterparts of Z sun and Xoebye- 
These are the only functions we need in the expansion for IT(P). 

In appendix C, we expand the necessary sum-integrals and integrals in 
powers of m/T. Substituting these expressions, keeping all the terms through 
order g 5 , the gap equation becomes 



2 1 2v/ 

m =2 9 t P 



1 m 2 

~p~2 p4 



+ I — 2 ~ l9 4 T 2 I Dchyc 

2 J p 2 + m 2 6 

p 



■yTf'iA-yTf^M. (193, 



2 y 7 p 2 + m 2 fpP A T J(p 2 + m 2 ) 2 ' 

p p 

The gap equation has poles in e and requires renormalization. Making the 
substitution g 2 — > Z g 2g 2 , the gap equation (193) can be written as 



A l9 2 \ 1 2 2 r l l 



+ T 



■ P 2 J p 2 + m 2 

v 

_I 4 T / 1 ^ J_ _ I 6 T f ^MP) 

r J p 2 + m 2 TpP i r J(p 2 + m 2 ) 2 ' 

p p 

(194) 

Using the expressions for the sum-integrals and integrals in appendix A and 
B and 
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a 2 y 
the gap equation (193) becomes finite: 



(195) 
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(1- log 2) a 2 



(196) 



Note that the running of the coupling in the three-loop ^-derivable approxi- 
mation agrees with the perturbative one (62) through order g A . 

Substituting the expansions of the sum-integrals and integrals into (183) 
and keeping all terms through order g 5 , the thermodynamic potential becomes 



+ m 2 ) - -T 
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p2 _|_ m 2 



~48 9 TpqrP 2 Q 2 R 2 S 2 2A 9 T P^ De¥ 
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(197) 



Using the expressions for the integrals and sum-integrals in the appendix, we 
obtain 



7T 2 T 4 
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(198) 



The divergences through order g 5 have been eliminated by using the gap equa- 
tion. In the calculations by Braaten and Petitgirard, they kept divergent terms 
through order g 7 in the gap equation and the thermodynamic potential. The 
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g(2TrT) 

Fig. 16. Two- and three-loop ^-derivable approximations to the pressure normalized 
to that of an ideal g clS clS cl function of g{2irT). Three-loop band is obtained by 
varying the renormalization scale, fx, by a factor of two around fj, = 2ttT. Note that 
using the running given by Eq. (173) at two-loops gives a result which is independent 
of fl. 

single poles in e of order g 6 could be eliminated by denning a renormalized 
coupling constant, but the double poles could not. The appearance of these 
double poles was identified with the truncation of interaction functional at 
three loops by Braaten and Petitgirard. However, it seems that it is possible 
to eliminate poles at all orders in the ^-derivable approach if the renormaliza- 
tion is performed properly [77] . Note that at the order presented here, namely 
g 5 , the difference between the two renormalization schemes does not enter since 
it only affects the treatment of terms of order g 6 and higher in the pressure. 

In Fig. 16, the two- and three- loop ^-derivable approximations to the 
pressure divided by that of an ideal gas are shown. Here we encounter the same 
problem as in SPT, namely that one cannot extend the solution to the gap 
equation past a critical value g* which we have indicated by a dot in the Figure. 
From the Figure it is evident that there is a significant improvement in the 
convergence compared to the weak-coupling expansion. Improved convergence 
properties were also found for the entropy and the Debye mass [80]. 

7.2 Gauge theories 

In this subsection, we review gauge theories within the ^-derivable ap- 
proach. This approach was first applied to gauge theories by Blaizot, Iancu 
and Rebhan [70,71,72] and by Peshier. [73] For simplicity, we restrict the dis- 
cussion to pure-glue QCD and take N c = 3 in which case the thermodynamic 
potential Q is given by 

Q[D, D gh ] = ^Tr log D' 1 - Tr log D^ 1 
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--Tr UD + Tr U gh D gh + P gh ] , 



(199) 



where the trace is over color and Lorentz indices. D(P) and P g h(P) are the ex- 
act gluon and ghost propagators, while II (P) and n gh (P) are the corresponding 
exact self-energies. $[P, D gh ] is the sum of all two-particle irreducible vacuum 
diagrams. 

The self-energy tensor U tlu (P) in Abelian gauge theories is transverse with 
respect to the four- momentum P = (u,p) and is also gauge-fixing indepen- 
dent. In nonabelian gauge theories, the self-energy tensor n M „(P) is generally 
not four-dimensionally transverse, but it can be decomposed into four different 
structure functions. [39] The axial gauges which are defined by n v A^ a = 0, 
where n M is a constant vector represent a class of gauges in which the self- 
energy tensor is transverse with respect to the four-momentum. However, an 
approximation does not necessarily respect transversality due to the fact that 
there are no vertex corrections. [72] In the temporal axial gauge § A = 0, 
this presents no problem since the self-energy tensor can be expressed exclu- 
sively in terms of and 11^. Moreover, in this gauge, the ghost self-energy 
n gh vanishes and the ghost field decouples. Thus the skeleton expansion of $ 
simplifies significantly. 




Fig. 17. One-loop QCD self-energy graphs. 

The gap equations that are obtained by varying the thermodynamic po- 
tential with respect to Py(P) and P^(P) are shown in Fig. 17: 

A T X = ui 2 - k 2 - n T , (200) 
A L 1 = fc 2 + n L , (201) 

where the transverse and longitudinal propagators are defined in Eqs. (148)- 
(149). The solutions to the nonlocal equations (200)-(201) are very difficult. 
However, to leading order in g and for soft momenta, the solutions are simply 
given by the hard-thermal loop self-energies (29)-(30) [72]. For hard transverse 
modes, Blaizot, Iancu and Rebhan solved the gap equations to next-to-leading 
order in g. The solution is 

n T = Tl {2) + 5T1 T , (202) 

5 At finite temperature, there are problems with this gauge related to the fact it 
is incompatible with the periodicity of the field. [98] There are other problems as 
well, see e.g. Ref. [111]. 
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where II^ is given by the one-loop self-energy diagrams with bare propagators, 
and SUt is given by the same diagrams with one line being soft and one line 
being hard. The calculation of the function SUt was performed in Ref. [72] 
where the explicit form can be found. 

Having found the approximate solutions to the gap equation, one can sub- 
stitute them into to expression (199) to obtain the pressure. In the imaginary 
time formalism, the only soft mode is the zero-frequency mode. In that case, 
the longitudinal self-energy II l reduces to a mass term m 2 D , while ITt van- 
ishes. For the nonzero Matsubara modes, it suffices to use bare propagators. 
The thermodynamic potential then reduces to 



7 



n = 4(d - l)f p log P 2 + 6(d - 1) V {$ p — 2 
+4T J log (p 2 + m 2 D ) . 



(203) 



Using the expressions for the sum-integrals and the integrals in the appendices, 
Eq. (203) reduces to 



8vr 2 T 4 r 15 a 15, 

1 1 I f3~~ \ A 

4 71 47T 3 
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(204) 



Using the weak-coupling expression for the Debye mass, it is easy to show that 
Eq. (204) agrees with the weak-coupling result for the pressure through order 
g 3 . The entropy, which follows from differentiation of (204) with respect to T 
becomes 



15 



45 I 4 7T 16tt 3 



{pm D y 



(205) 



The expression for the entropy (205) also agrees with the weak-coupling ex- 
pansion result to order g 3 . 

Again Blaizot, Iancu and Rebhan [72] exploited the fact that the two- loop 
expression for the entropy reduces to the one-loop expression. In the case of 
gauge theories, one can write [72] 



<S — Sl + <St , 
where 

S L = -8 



-16 



f du 


r dn{uj) 


1 2^ J 


dT 



du f dn{uj) 



J 2^ J 



dT 



Imlog (k 2 + II L ) + Im II l Re D l 
Im log (-w + k 2 + II T ) - Im I1 T Re D T 



(206) 



(207) 



(208) 
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The Stefan-Boltzmann result for the entropy is obtained by using bare 
transverse and longitudinal propagators in Eqs. (207)-(208). The order-g 2 
correction to the entropy is obtained by using the leading-order solution to 
the gap equation in Eq. (208). Note that the longitudinal modes do not 
contribute at this order in g. The order-g 3 correction receives contribution 
from both the soft and hard modes. The soft contribution can obtained from 
Eqs. (207)-(208) by using the hard-thermal-loop self-energies and expanding 
the Bose-Einstein distribution function, n(u>) ~ T/uj. For example, the con- 
tribution from soft longitudinal soft gluons is 



The hard contributions are found by evaluating Eq. (208) using the next- 
to-leading order solution to the gap equation for transverse modes 5Ut(ou, k) 
on the light cone, i.e. uj = k: 



Blaizot, Iancu and Rebhan define the HTL entropy by using the hard- 
thermal-loop self energies in Eqs. (207)-(208) for both soft and hard momenta. 
While this is not the leading-order solution to the gap equation for the hard 
modes, it nevertheless reproduces the order-g 2 term in the entropy (since this 
is obtained by using the one-loop self-energies with bare propagators in the 
expression for the entropy and is exclusively given by the hard modes). The 
full numerical evaluation of the HTL entropy is of course also nonperturbative 
since its expansion contains all orders in the coupling g. 

The ultraviolet divergences that appear in the approximate solutions to the 
one-loop gap equations can be eliminated by wave function renormalization at 
zero temperature. However, they drop out of the expression for the entropy. 
Thus consideration of the entropy gives no renormalization group equation 
describing the running of the coupling constant. In their numerical treatment, 
Blaizot, Iancu and Rebhan took the standard two-loop running of the coupling 
constant. The HTL resummation in Refs. [70,71,72] is specific to the entropy. 
In order to obtain the pressure, one must integrate the entropy with respect 
to the temperature. This introduces an unknown integration constant which 
the authors fixed by requiring that the pressure vanish at T c . 

Peshier [73] considered another approximation to the two-loop <l>-derivable 
approximation for gauge theories motivated by scalar field theory. The two- 
loop ^-derivable approximation is equivalent to making the substitution 




(209) 



k 




(210) 



k 




(211) 
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where m is the solution to the one-loop gap equation (168). If one replaces 
the mass m by the hard-thermal-loop approximation to the self-energy, one 
reproduces the pressure correctly to order g 2 . The order-g 3 contribution is, 
however, overestimated by a factor of four. Peshier suggested to approximate 
the two- loop approximation to the functional <&[D] in QCD by the analog to 
Eq. (211): 

$[£>*] = ^TiD*U* , (212) 

where the * indicates that the self-energies and propagators are taken in the 
HTL approximation. By construction, this approximation to the two-loop $- 
derivable thermodynamic potential reproduces the order-g 2 correction to the 
pressure. The order-g 3 is again incorrect by a factor of four. The entropy of 
Peshier is calculated by differentiating (212) with respect to the temperature. 
The resulting entropy does not coincide with that of Blaizot, Iancu and Reb- 
han. 

Finally, we would like to discuss the strategy developed by Braaten and 
Petitgirard to solve the n-loop ^-derivable approximation in QCD. We know 
that the static infrared limit of n 00 (u;, p) is the order-g 2 result for the Debye 
mass m 2 D , while all other components vanish. For soft momenta, we therefore 
expand 

U 00 (0,p)=m 2 D + g i [a 4 ' ~ 2 (p) + -] , (213) 

where mo is the Debye mass is defined by the pole position of the static 
propagator: 

m 2 D = U 00 (u n = 0,p)\ p=tmD . (214) 
For hard momenta, we write the expansion of n 00 in the form 

IMP) =g 2 [if ' °( P ) + rC(p) + ...] . (215) 

In order to solve the gap equation (214) beyond leading order, we must 
incorporate the self-energies in the propagators. In a next-to-leading-order 
calculation of the pressure, it suffices to use bare propagators for hard loop 
momenta: 



5 P P 

p 2 ~T~ \ L W p 4 



(216) 



where £ is the gauge parameter in the class of covariant gauges. For soft 
momenta, we include the leading term in the expansion for n o- We must 
then use the dressed propagator: 
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For soft momenta, the other components of the propagator are still given 
by (216). The gap equation then becomes: 



ml = 3(d-l)yf K ±- 2 +3g>Tj 



d-2 + 2(m 2 D -p 2 ) 



k 2 k 2 [(p + k) 2 + m 2 D ] 



k 2 + m 2 D vs J ^ u > k 4 [(p + k) 2 + m 2 D 



P 2 =-m 2 D 



. (218) 



The first term arises when all internal lines are hard and is independent of 
the gauge parameter £. The remaining terms arise when all lines are soft. The 
case where one line is soft and one line is hard precisely corresponds to 5Ut 
in (202) and is neglected in a NLO calculation. The contribution from the soft 
sector was first calculated by Rebhan when he calculated the order g 3 -term 
in the weak-coupling expansion for the Debye mass using resummed pertur- 
bation theory [113]. When this contribution is evaluated on shell, the gauge 
dependent term vanishes algebraically before we have integrated over k. How- 
ever, the integral is linearly infrared divergent. Introducing an infrared cutoff, 
performing the integral and taking the limit p 2 — > —m 2 D before removing the 
cutoff, it can be shown that this term drops out [112,113,114]. 

The third term is logarithmically infrared divergent and this singularity is 
caused by a massless static gluon propagator. This is another example of the 
breakdown of perturbation theory in the static magnetic sector of QCD and 
prevents a perturbative definition of the Debye mass beyond leading order in 
g 2 . A magnetic mass m m of order g 2 T is expected to be generated nonpertur- 
batively in QCD. Rebhan [113] suggested to obtain an estimate of the NLO 
correction by incorporating m m by the replacement 1/p 2 — > l/(p 2 + ttl 2 ^ in 
the transverse part of the propagator. The integrals can then be evaluated 
analytically and the gap equation (218) then becomes 



2 2^2 . 3g 2 m D T 
m D = gT + 



1 _ m ^ \ lo 2m D + m m _ 1 _ m m 
Am 2 D ) m m 2 2m D 



, (219) 



where the first term comes from hard loop momenta and the remaining terms 
from the soft sector. The pressure is still given by Eq. (204), where the mass 
mo is the solution to (219) instead of the weak-coupling expression that was 
used in Ref. [130]. By introducing a magnetic mass, we have included nonper- 
turbative physics in an ad hoc manner. One might abandon the NLO solution 
to the gap equation altogether and simply use the leading-order result. Our 
approximate solution to the two-loop ^-derivable approach then reduces to 
the leading-order solution of Blaizot, Iancu and Rebhan [72]. 
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8 Dimensionally Reduced SPT 



In Sees. 5-6, we have seen that SPT and HTLPT are ways of reorganiz- 
ing perturbation theory so that the convergence properties of the perturba- 
tive approximants is dramatically improved. However, we have also seen that 
the two-loop result for the HTL-improved pressure in pure-glue QCD fails to 
match onto that of lattice gauge theory at temperatures for which they are 
available. It has been suggested that the disagreement is due to the fact that 
the HTL approximation breaks down for hard momenta. One should therefore 
use the HTL approximation only for soft momenta and treat hard momenta 
in a different manner. Such a strategy has recently been proposed by Blaizot, 
Iancu and Rebhan [130], where one uses strict perturbation for the hard modes 
and applies screened perturbation theory only to the dimensionally reduced 
theory for the soft modes. 

This approach is called dimensionally reduced screened perturbation the- 
ory (DRSPT). In the case of gauge theories, this approach is gauge invariant 
order by order in the loop expansion. Of course, for nonabelian gauge theo- 
ries, it breaks down at four loops due to infrared divergences in the magnetic 
sector. This is reflected in Eq. (125) where the constant a can only be deter- 
mined nonperturbatively. Blaizot, Iancu, and Rebhan [130] applied DRSPT 
to EQCD and calculated the pressure through four loops up to this unknown 
constant. It seems that the agreement with available lattice data is improved 
significantly by treating the hard modes perturbatively but not insisting on 
expanding the soft contributions in powers of the coupling, instead keeping 
the full (^-dependence of these contributions. 

In this section, we apply DRSPT to a dimensionally reduced scalar field 
theory. The effective Lagrangian is then written as 

C= i(V0) 2 + \{m\ + 5m 2 )0 2 + ^ - Um^ 2 , (220) 

where is the mass parameter determined perturbatively in the dimensional 
reduction step and 5m is to be determined variationally: 

— ^/ = °- (221) 



d(5rri 

In the four-loop calculation below, we also need the coupling constant A to 
next-to-leading order in g 4 : 

The diagrams through four loops are shown in Fig. 18. 
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1- loop: |Q 

2- loop: iQQ-lQ 

3- .00P: -iOD-AOOO+iCO-iO 

4- .OOP: ,@ + ,(^0 +ic&+ xOC)CO 

-iOOO-iOOO-AdD 
+KXMOCHG 

Fig. 18. Diagrams which contribute up to four-loop order in DRSPT. A boldfaced 
x indicates an insertion of 5m 2 . 

8.1 One-loop order 

The one-loop contribution to the free energy is 

f = l -J\ og (p 2 + 



m 2 



p 

1 ,3 



127T 



m 6 . (223) 



with m 2 = m\ + 5m 2 . At this order the only solution to the variational 
condition dfo/d(5m 2 ) = is the trivial solution 5m 2 = or m 2 = m 2 E . 

8.2 Two-loop order 

The two-loop contribution to the free energy is 

2 

5m 2 



fl 8 A (/p 2 + m 2 ) 2 



p 



pi _|_ m 2 



i Tfi 

= -am 2 T + —5m 2 , (224) 

8 07T 

where we have used the integrals contained in appendix B and again m 2 = 
m\ + 5m 2 . Combining the one- and two-loop contributions and requiring that 
their sum satisfy the variational condition, gives the two-loop gap equation 

m 2 - m \ = -2namT . (225) 

The two-loop DRSPT approximation to the free energy is given by the 
sum of (87) truncated at order a, (223) and (224). Using the leading-order 
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solution to the gap equation (225), it is easy to show that this approximation 
agrees with the weak-coupling expansion through order g 3 . 



8.3 Three-loop order 



The three-loop contribution to the free energy is 

/2 = _1 A 2 f _± 1 1 1 

48 J p 2 + m 2 q 2 + m 2 r 2 + m 2 (p + q + r) 2 + m 2 

pgr 



6 A ' (/p2 +m2 ) /( 2 + m 2)2 + 4 A5m2 /(p 2 + 
\p J q pq 



16 



1 



(p 2 + m 2 ) 2 dm 2 



where the mass counterterm is given by 



m 2 ) 2 q 2 + m 2 
(226) 



A 2 m 2 



A 



3 y 167T J e 

Using the integrals listed in Appendix B, Eq. (226) reduces to 
1 



(227) 



f 2 = -na 2 mT 2 
3 



A 9 
log 1 log 2 

& 2m 8 . 



-aT5m 2 - ' ' 2 - 



3271171 



(5m 2 Y . (228) 



Combining this result with the one- and two-loop contributions (223)-(224), 
and imposing the variational constraint gives the three-loop gap equation 



(m 2 — 



rn 



a 2 m 2 T 2 



l0g 2^- l0g2 + ^ 



(229) 



The three-loop DRSPT approximation to the free energy is given by the sum 
of (87), (223), (224). and (228). Upon substituting the expression for the mass 
parameter (83) and using (229), the A-dependence in the three- loop DRSPT 
approximation to the free energy cancels to next-to-leading order in g. This 
approximation agrees with the weak-coupling expansion to order g A . 

8.4 Four-loop order 

The four-loop contribution to the free energy is 

A 3 r 1 1 1 1 1 



/3= A^ f ill 1 

48 J p 2 + m 2 q 2 + m 2 r 2 + m 2 s 2 + 

pqrs 



m 2 (p + q + r) 2 + m 2 
1 



X 



(p + r + s) 2 + m 2 



63 



0.98 



| 0.96 



0.94 



0.92 



0.90 




l_\ DRSPT Two Loop 
L_J DRSPT Three Loop 
1 DRSPT Four Loop 



2 

g(2nT) 



0.98 



S 0.96 



0.94 



0.92 



0.90 



l_\ Perturbative 
L."J DRSPT Three Loop 
| SPT Three Loop 



2 

g(2rcT) 



Fig. 19. (a) Two-, three-, and four-loop DRSPT results for the pressure as a func- 
tion of g(2irT). (b) Three-loop DRSPT and SPT results are compared with the 
perturbative result accurate to order g 5 . 
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+ |^A 2 m 2 + 5/ 3 , 



(230) 



where A 2 m 2 is given in Eq. (227) and 5/3 is a vacuum counterterm given by 



5/3 



7T 4 A 3 1 

96 (4tt) 6 e 



(231) 



The logarithmic ultraviolet divergence that is cancelled by the vacuum coun- 
terterm (231) is arising from the first four-loop diagram in Fig. 18. It is the 
same divergence that we discussed in Sec. 4.1.2 that leads to the running of 
/(A) at four loops. Thus depends explicitly on A at order A 3 and this de- 
pendence is cancelled only by adding the hard contribution to the free energy 
through order g 6 . 

Using the integrals in Appendix B, we obtain the four-loop contribution 
to the free energy 
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h = 



7T 2 (7T 2 -4) 



a 3 T 3 1()g 0.845213 



12 



2m 



6m 



7T 




192vrm 3 



1 



(5m 2 ) 3 . (232) 



Note that the contribution proportional to a (5m 2 ) vanishes. Combining this 
with the lower order results and imposing the variational condition gives the 
four-loop gap equation 



Note that Eq. (233) contains a contribution from the two-loop contribution 
(224) evaluated using the NLO expression for A given by Eq. (222). The soft 
part of the four-loop DRSPT approximation to the free energy is given by the 
sum of (223), (224), (228) and (232), again with (224) evaluated using the 
NLO expression for A given by Eq. (222). 

In Fig. 19 we have plotted the two-, three-, and four-loop DRSPT results 
for the pressure as a function of g(2-n;T). The fact that the four-loop DRSPT 
approximation deviates from the two and three-loop DRSPT approximations 
is related to the fact that we have taken the hard contribution to the free 
energy of order a 3 to be identically zero which is highly improbable; however, 
absent of a calculation of this coefficient we prefer to not play games. Un- 
fortunately, the calculation of the contributing four-loop diagrams is a highly 
nontrivial task. Note that to obtain the three-loop approximation we used the 
gap equation obtained at four-loops which is independent of this unknown 
constant. We use this gap equation because the three-loop gap equation does 
not have a solution for g ^ 0. ^ In Fig. 19(b) we compare the three-loop 
DRSPT result with the g 5 perturbative result and the three-loop SPT result. 
As can be seen from this Figure, although DRSPT improves the situation in 
terms of the convergence of successive approximations, the result is still very 
sensitive to the renormalization scale having a larger variation than even the 
perturbative result. In the next section we will demonstrate an alternative 
method for reorganizing the dimensionally-reduced calculation. 

9 Dimensionally Reduced ^-derivable Approximation 

In this section, we apply the ^-derivable approach to the dimensionally 
reduced (DR) scalar field theory defined by (69). The total free energy is given 
by the sum of the hard and the soft contributions. The hard contribution 

' This is similar to the behavior found in the three-loop SPT gap equation [88] and 
the three-loop gap equation in DRSPT applied to QCD. [130] 




- — a 2 m 2 T 2 (m 2 - m 2 F ) log log 2 + - 

3 V EJ [ 6 2m 6 8 




65 



is given by strict perturbation theory and is an expansion in powers of a. 
Through order a 2 , the expression is given by (87). The soft contribution is 
given by calculating the free energy for the effective theory (69) in the three 
dimensions using the ^-derivable approach. We will use the strategy developed 
by Braaten and Petitgirard in Ref. [80] and explained in Sec. 7. 



9.1 Two-loop DR ^-derivable approximation 



The thermodynamic potential il 2 in two-loop DR ^-derivable approxima- 
tion is given by 

i r n(p) 



fi 2 = i|log(p 2 + m| + n(p)) 



+ l x (j ^^TW) 1 • (234) 

Again the self-energy Tl(p) is momentum independent. If we write m 2 = 
m 2 E + IT, the gap equation that follows from varying (234) is 



2 J p 2 + m 2 E + Yi(p) 

2 



1 r 1 

7- A / - . (235) 

2 J v 2 + m 2 y ' 



2 2 , 

m =m p + 

" //- • /// 

p 



Both the thermodynamic potential and the gap equations are finite with di- 
mensional regularization. Using (B.2) and (B.3) in appendix B, we obtain 

n 2 = -^--±m 2 Ta, (236) 
m 2 = m E - 2namT . (237) 

Note that the gap equation and the thermodynamic potential are the same 
as those obtained in the two-loop DRSPT approximation. The two-loop DR $- 
derivable approximation then agrees with the weak-coupling expansion through 
order g 3 . 



9.2 Three-loop DR $- derivable approximation 



The three-loop DR ^-derivable thermodynamic potential f2 3 is given by 

n(p) 



^ = \j\og{p 2 + m 2 E + Ii(p))- l -j 



v 



p 2 + m 2 E + n(p) 



8 \J p 2 + m 2 E + U(p)j 48 J qr p 2 + m\ + n(p) q 2 + m\ + U(q) 
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X r 2 + m| + n(r) s 2 + m\ + II(s) + 2 A ' m ^ / p 2 + m| + II(p) ' (238) 
where A 2 m E is the mass counterterm and s = — (p + q + r). 

A 2 m 2 = — ^ — — . 239 

This mass counterterm is identical to the one used in the dimensional 
reduction step and in DRSPT. The gap equation obtained by varying (238) 
with respect to II (p) is 



q 2 + m\ + U(q) 



6 A ' / q 2 + m 2 E + 11(g) r 2 + m| + II(r) s 2 + m| + U(s) + A2 "^ ' (240) 



qr 



fl 3 can be simplified by substituting the gap equation (240) into (239): 

n(p) 



fi 3 = ljlog(p 2 + m| + II(p))-ij 



4 7 » 2 + m F + n(p) 



48 ^ p 2 + m| + II(p) g 2 + m| + 11(g) 

hm + I A2m| / , 2 + m| + n( P ) • (241) 



1 

X 



r 2 + m| + n(r)s 2 + m| + Il(s) 4 ^ 7 p 2 + m| + II(p) 

In the three-loop DR ^-derivable approximation, the gap equation is no longer 
independent of the external momentum and we shall carry out the same mass 
expansion as we did in Sec. 7. The Debye mass now satisfies 

m 2 = m 2 E + \\f I - \ A 2 / Debyc + A 2 m| . (242) 

2 J Q + n^E + n(g) 6 

The variational equation (240) then reduces to 

n(p) = m 2 - m| + ^A 2 [J Dcbyc - J sun ] . (243) 

We distinguish between hard and soft loop momenta and the expansions for 
the self-energies are, respectively, 

n(p) = m 2 - m 2 E + A 2 [n 4 , (p) + 114,1 (P) + •••] + -, [P hard] , (244) 
n(p) =m 2 - m\ + A 2 [<t 4 _ 2 (p) + (T4,o (p) + ...] + ..., [p soft] . (245) 
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It should be noted that the expansion (244) for hard momenta is never needed. 
The reason is that if one momentum p in a loop integral is hard, one can expand 
the self-energy in powers of Yl{j>)/p 2 . The resulting integrals have no scale and 
vanish in dimensional regularization. H 

By inserting the expansion for a(p) into the gap equation, we can deter- 
mine the functions <7 n k (p) by matching coefficients of order A n on both sides. 
The only function we need to carry out calculations to order g 5 , is a^^ip) 
which is given in Eq. (192). The gap equation reduces to 



2 2 1 x f 1 1,2, 
m ~ m E = o A / — t 2 ~ r A /o ebye 

2 J v + m z 6 



V 



p 

Using the expressions for the integrals in appendix B and for the mass coun- 
terterm (239), the renormalized gap equation reduces to 



m 2 — m 2 F = —2namT n 2 a 2 T 2 

E 3 



( log 2^ + |- 21 ° s2 ) 



T 

+vra— (1 - log 2) 
m 



(247) 



After substituting the expression (83) for the mass parameter, the A-dependence 
in (247) cancels to next-to-leading order in g. 

Again, using the expansions of the integrals in appendix C, the dimension- 
ally reduced three- loop ^-derivable approximation, f2 3 , to the thermodynamic 
potential through order g 5 reduces to 

~ 1 f, / 9 9n 1 fm 2 —m% 1 , 9t f 1 



V V 



_1 A 2 f il l 1 

48 J p 2 + m 2 q 2 + m 2 r 2 + m 2 (p + q + r) 2 + m 2 

pgr v ' 

+7 AW / + , / • ( 248 ) 

4 J (p 2 + m 2 ) 2 4 J p 2 + m 2 

p p 

Using the expressions for the integrals in appendix B and for the mass coun- 
terterm (239), we obtain the soft contribution to the renormalized three-loop 
free energy 



m 3 m 



i2, 16, ( m2 - m °) + r mTV ( log 4 + 2 - 5 log 2 ) • (249) 



" This is in contrast to the ^-derivable approximation in 3+1 dimensions, where 
the nonzero Matsubara frequencies set the scale. 
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g(27rT) 

Fig. 20. DR <£-derivable approximations to the pressure divided by that of an ideal 
gas. For comparison we have shown the three-loop SPT result for the pressure. 

The complete three-loop DR ^-derivable approximation to the free energy 
is given by adding (87) and (249). Upon substituting the expression for the 
mass parameter (83) and using the fact that m 2 is independent of A to next- 
to-leading order, it is easy to show that A-dependence in the soft contribu- 
tion (249) cancels to next-to-leading order in g. Using Eq. (247), one can show 
that the complete three-loop DR ^-derivable approximation to the free energy 
agrees with the weak-coupling expansion through order g 5 . 

In Fig. 20, we show the two- and three-loop DR ^-derivable approxima- 
tions to the thermodynamic pressure. The bands are obtained by varying the 
renormalization scale by a factor of two around the central value \i = 2nT. 
From the Figure, we see that DRPHI also has better convergence than the 
perturbative expansion (61). Note that for \x = 4-irT the DRPHI three-loop 
gap equation has no solution beyond g ~ 3.2 which we have indicated by 
not continuing the topmost DRPHI line in Fig. 20. For comparison, we have 
also shown the three-loop SPT approximation to the pressure. The agree- 
ment between the two approaches is very good and the final scale variation 
of the DRPHI result is smaller than the SPT result which is most likely due 
to the fact that with three-loop SPT we are forced to use a two-loop gap 
equation owing to the sickness of the three-loop variational gap equation. The 
dimensionally reduced ^-derivable approach therefore offers an economical al- 
ternative to the full treatment in four dimensions. However, as in the case of 
DRSPT, it is limited to calculating static quantities. 

10 Conclusions 

In this paper we have reviewed some of the progress that has been made 
concerning our understanding of ultrarelativistic hot field theories over the 
past 10-15 years concentrating on the systematic computation of thermo- 
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dynamic observables. We have presented results using multiple methods of 
computation and perturbative reorganization. In particular we have discussed 
standard perturbative techniques, dimensional reduction, effective field theory 
methods, screened perturbation theory (SPT), hard-thermal-loop perturba- 
tion theory (HTLPT), the ^-derivable approximation, dimensionally reduced 
(DR) SPT, and the DR ^-derivable approximation. 

We began by showing that the perturbative expansions of the free energy 
for both simple scalar field theories and QCD have a very small, if not van- 
ishing, radius of convergence. In the case of QCD this makes the resulting ex- 
pansions useless at phenomenologically relevant temperatures where a s ~ 0.2 
(g s ~ 2). We then showed that dimensional reduction and effective field theory 
methods can provide an efficacious way of organizing weak-coupling calcula- 
tions. By combining them with nonperturbative methods such as lattice gauge 
theory, the long-standing problem of calculating the order-a^ contribution to 
the QCD pressure can be solved. 

Another advantage of using dimensional reduction and effective field the- 
ory methods is that the explicit separation of the hard and soft scales al- 
lows for a different treatment of the physics at these scales. We showed, for 
example, that if the soft-sector is treated non-perturbatively, the successive 
approximations to the thermodynamic functions are much better behaved, 
giving reasonably constrained results for the pressure in scalar field theories 
to g ~ 4. We showed that two reorganizations of the soft-sector computation, 
the DRSPT and DR ^-derivable approach, both gave better convergence of 
the successive approximations to the pressure at large coupling. 

While the results are impressive, the methods are limited because they 
can only be applied to the computation of static quantities such as the pres- 
sure. It would be preferable to have a reorganization of the soft-sector physics 
which could be applied to the calculation of real-time quantities as well. This 
desire has motivated the development of the ^-derivable approach and the 
SPT and HTLPT reorganizations. The ^-derivable approach has only been 
taken to three- loops for a scalar field theory where it improves covergence 
and scale variation of the approximations to the pressure. It would be nice to 
systematically extend this approach also to three loops for QCD. It has been 
demonstrated that when applied to the calculation of the pressure in QCD 
the HTLPT reorganization yields impressive results when one considers the 
convergence of the successive approximations. Additionally, results are explic- 
itly gauge-invariant within HTLPT as opposed to ^-derivable approximations 
where gauge invariance can only be guaranteed to a fixed order in the coupling 
constant. 

Despite these theoretical benefits, when compared to existing lattice data 
for the pressure in the range T c < T < 5T C , the NLO HTLPT prediction 
for the pressure seems to underestimate the correction to ideal gas behavior. 
The failure of the method in this region could come from the fact that within 
HTLPT the requirement of gauge invariance causes a simultaneous modifica- 
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tion of both soft and hard modes. Perhaps a method that makes an explicit 
separation between these two scales is more appropriate. Another possible ex- 
planation is that HTLPT, as implemented, does not properly take into account 
the proper scaling relations for quantities near the phase transition which is 
governed (approximately) by the center group of SU(3), Z(3). However, de- 
ciding which of these two possibilities is the culprit is not possible at this 
stage. Results from the approximately self-consistent ^-derivable approach of 
Blaizot, Iancu, and Rebhan seem to point to the first possibility, but are by 
no means conclusive. The only thing that can resolve this is further study. We 
also mention in closing that given the success of the DR ^-derivable approx- 
imation discussed in Sec. 9 it would be interesting to apply this method also 
to the QCD pressure. 
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A Sum-integrals 

In the imaginary-time formalism for thermal field theory, the 4-momentum 
P = (PoiP)is Euclidean with P 2 = Pq + p 2 . The Euclidean energy po has 
discrete values: P — 2n7rT for bosons and P = (2n + l)nT for fermions, 
where n is an integer. Loop diagrams involve sums over P and integrals over 
p. With dimensional regularization, the integral is generalized to d = 3 — 2e 
spatial dimensions. We define the dimensionally regularized sum-integral by 

' eVV r E f£&- (a.i) 



1 

{P} 



47r / Po^T J ( 27r 

e W <r -'l> 

P = (2n+l)wT 



| .. ! T E J ( 27 r)S-2e ' ( A - 2 ) 

^ ) Pn = C2n+lWT J \ Z7r ) 



where d = 3 — 2e is the dimension of space and /i is an arbitrary momentum 
scale. The factor (e 7 /47r) e is introduced so that, after minimal subtraction of 
the poles in e due to ultraviolet divergences, fj, coincides with the renormal- 
ization scale of the MS renormalization scheme. 

A.I One-loop sum-integrals 

We need several simple massless one-loop sum-integrals 
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logP 2 = --T 4 + 0(e), 
p 45 

1 / H \ 2t 1 

pP 2 ~" \J^f) 12 



C(-i) 



2e r 



p(P 2 ) 2 (4tt) 2 V4ttT 
r 1 T 2 



+ 2 7 + 0(e) 



(A.3) 
(A.4) 

(A.5) 

(A.6) 



The calculation of these sum-integrals is standard. We also need some massive 
one-loop sum-integrals 



rn 



2c 



1 



(47r) 2 \m. 
1 / n n 2e 



pP 2 + m 2 (An) 2 \m 
1 1 / /i n 2e 



p(P 2 + m 2 ) 2 (4tt) 2 Vm 



e(l-e)(2-e) 

e(l-e) 
e 7e r(l + e) 



+ J 2 



, (A.7) 
(A.8) 
(A.9) 



The thermal terms can be expressed as integrals involving the Bose-Einstein 
distribution function: 



J n ((3m) 



' 12 j (3 4 - 2n m^Jdk * 



1 



r(§-n-c 







(k 2 + m 2 ) 1/2 eP^+ m2 ) 



1/2 



These integrals satisfy the recursion relation 

xJ' n (x) = 2tJ n (x) - 2x 2 J n+l (x) . 
In the limit (3m — ► 0, these integrals reduce to 



(A.10) 



(A.ll) 



Jn 



16vr 4 4tt 2 yj 8tt 3 
— — (pm) + — (pm) 



Ji 
J 2 



45 



47T 2 

3 

2tt 



4ir/5m - 2 hog ^ - i + 7 J (/3m) 2 , 



(A.12) 
(A.13) 
(A.14) 



A. 2 Two-loop sum-integrals 



We need a single massless two-loop sum-integral which happens to vanish: 
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¥ pq p 2 q 2 (p + q) 2 °' 

A. 3 Three-loop sum-integrals 

We need a single massless three-loop sum-integral: 
1 T 4 



(A.15) 



pqrP 2 Q 2 R 2 (P + Q + R) 2 24(4tt) 2 



(A.16) 



B Three-dimensional integrals 

Dimensional regularization can be used to regularize both the ultraviolet 
divergences and infrared divergences in 3-dimensional integrals over momenta. 
The spatial dimension is generalized to d = 3 — 2e dimensions. Integrals are 
evaluated at a value of d for which they converge and then analytically con- 
tinued to d — 3. We use the integration measure 



47T 



d 3 ~ 



2c 



P 



(27T) 



3-2e 



(B.l) 



We require several integrals that do not involve the Bose-Einstein distribution 
function. The momentum scale in these integrals is set by the mass parameter 
m or mo- 

B.l One-loop integrals 



The one-loop integrals are 

,3 



log [p 2 + m 2 ) 
1 



/ 

p 

/ 



m 
6n 

m 



p 2 + m 2 47r 



(p 2 + m 2 ) 8nm 



(p 2 + m 2 f 327rm 3 



0-4 -2 (p) 



l+(21og^ + 2) e + 0( e 2 ) 



' !21og^) e + 0(e 2 ) 



I i21og^-2;. 0(r^ 



T 2 



(p2 + m 2)2 i2(47r) 3 m 



;i-log2) . 



(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 
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The calculation of (B.2)-(B.4) are standard. The integral in (B.6) was calcu- 
lated by Braaten and Petitgirard [80]. 

B.2 Two-loop integrals 



The only two-loop integral required is 

1 



'Debye 



/ 



qr 



4(4tt) 2 

B.3 Three-loop integrals 



(q 2 + m 2 )(r 2 + m 2 )((p + q + r) 2 + m 2 ) 



- + 41og-^ + 6-81og2 + (9(e) 
e 2m 



p=tm 



The three-loop integrals required are 

/■111 1 
J p 2 + m 2 q 2 + m 2 r 2 + m 2 (p + q + : 



pgr 



m 



1 



1 



(4tt) 3 

1 



r) 2 + m 2 



- + 6 log + 8 - 4 log 2 + (9(e) 
e 2m 



pgr 



^2 _|_ m 2 g2 _|_ m 2 r 2 _|_ m 2 _|_ q _|_ r ^2 _|_ m 2ji 



8m(47r) 3 
B.4 Four-loop integrals 



I + 6 log-^ + 2-41og2 + C(e) 
e 2m 



The only four-loop integral needed is 
1111 



1 



pqrs 



p 2 + m 2 q 2 + m 2 r 2 + m 2 s 2 + m 2 (p + q + r) 2 + m 2 

1 



(p + r + s) 2 + m 2 



7T 



32(4tt) 



I + Slog A + 2 + 41og2-^C(3) + 0(e) 



2m 



7T^ 



This was calculated in Ref. [107]. 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



C m/T expansions of sum-integrals 

We next consider the contribution to the free energy F 2 h from the bas- 
ketball diagram. There are two momentum scales in that the sum-integral; 
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the hard scale T and the soft scale gT. The hard scale includes all nonzero 
Matsubara frequencies for all momenta p and also the region p of order T 
with n — 0. The soft scale involves the region p of order gT with n — 0. Since 
there are three sum-integrals, all momenta can be hard, two momenta can 
be hard and one soft, two momenta can be soft and one hard, or all three 
momenta can be soft. These regions are denoted by (hhh), (hhs), (hss), and 
(sss), respectively. If a momentum is hard, one can expand in powers of m/T. 
If a momentum is soft, one must keep the mass in the propagator. We begin 
by writing 



J~2\> — — ^g^ball , 

where Xbaii is the basketball sum-integral: 

1 1 1 



(CI) 



-ball : 



1 



' pqrP 2 + m 2 Q 2 + m 2 R 2 + m 2 S 2 + m 2 ' 
where S — — (P + Q + R). The various contributions to Z ba ii are 

1 



(C.2) 



r (hhh) 

'ball 



-i 



q-{hhs) 
X ball 



j-(hss) 
'ball 



AT 



pqrP 2 Q 2 R 2 (P + Q + R) 2 ' 

1 



I p 2 + m 2 



AT 1 



1 



qrQ 2 R 2 {Q + R) 2 ' 
1 _f' 1 



pq 



T (SSS) _ rp3 

X ball — J 



p 2 + m 2 q 2 + m 2 
1 1 



R R 4 ' 
1 



pqi 



p 2 + m? q 2 + m 2 r 2 + m 2 (p + q + r) 2 + m 2 



(C.3) 
(C.4) 

(C.5) 

(C.6) 



Using Eqs. (A. 15), (B.3), and the fact that m ~ gT, we obtain Eq. (54). 

We need the m/T expansions of various one- loop sum-integrals appearing 
in the two and three-loop ^-derivable approximation to the free energy. We 
need the expansions for both soft and hard momenta. Inserting the expan- 
sions (188) and (189), we obtain for the one-loop sum-integrals 



^log (p 2 + n(P))=^ 



m 1 m 
T 2 ~ 2 P 1 

Kf n 4)0 (P) + n 4)1 (p) 



logP 2 + 

+g 



p 2 



^ log (P 2 + n(P)) = T J log {f + m 2 ) + g 4 T J °- 



- cr 4 _2 (p) 



+ m 1 



,w n(p) 



fp P 2 + n(P) 





m 


41 

m 






~ pi 



(C.7) 
(C.8) 



+ ^_ n 4 , o( P) + n 4 , l( P) (C9) 
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w n(P) 



p p 2 + n(P) 
1 



m 



p P 2 + n(P) ./ /> 

w 

p p 2 4 



^2 _|_ m 2 

1 m 2 

p2 p4 



f P P 2 + U(P) T Ip 2 + m 2 9T I 



(p 2 + m 2 ) 2 



04,-2 (J?) 

(p 2 + m 2 y 



(CIO) 

(C.11) 
(C.12) 



Similarly, we need the m/T expansions of some two-loop sum-integrals 



f (hh) 
-Debye 

f (hs) 
'Debye 



1 



QR Q 2 R 2 (Q + R) 2 



3¥] -tT 



p P 4 J p 2 + m 2 



Olg 



T (SS) _rp2j 

X Debye ~ 1 1] 



Debye 



(C.13) 
(C.14) 

(C.15) 



Finally, we need the m/T expansion of a single three- loop sum-integral. Through 
order g 5 we have 



1 



1 



1 



Ahhh) I 

Tpqr p 2 + n(P) q 2 + n(Q) p 2 + n(p) s 2 + u{S) 

PQR p 2 Q 2 R 2 s 2 +0 ( g2 ) ' 

1111 



f(sss) 

Tpqrp 2 + n(P) q 2 + u{Q) r 2 + n(p) s 2 + u(S) = 
[ 1 1 1 1 o ( 2 ) 

J p 2 + m 2 q 2 + m 2 r 2 + m 2 (p + q + r) 2 + m 2 W ' 

pgr 

The (hhs) and (/iss) terms first start to contribute at order g 2 . 



(C.16) 



(C.17) 
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